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Abstract 

Here we review and extend central limit theorems for highly chaotic but deterministic semi- 
dynamical discrete time systems. We then apply these results show how Brownian motion-like 
results are recovered, and how an Ornstein-Uhlenbcck process results within a totally determin- 
istic framework. These results illustrate that the contamination of experimental data by "noise" 
may, under certain circumstances, be alternately interpreted as the signature of an underlying 
chaotic process. 
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1 Introduction 

Almost anyone who has ever looked through a microscope at a drop of water has been intrigued by 
the seemingly erratic and unpredictable movement of small particles suspended in the water, e.g. 
dust or pollen particles. This phenomena, noticed shortly after the invention of the microscope by 
many individuals, now carries the name of "Brownian motion" after the English botanis t Robert 
Brown who wrote about his observations in 1828. Almost three-quarters of a century later, EinsteinI 
(190i) gave a theoretical (and essentially molecular) treatment of this macroscopic motion that 
predicted the phenomenology of Brownian motion. (A very nic e Engl i sh tra nslation of this, and 
other, work of Einstein's on Brownian motion can be found in iFiirthI (jl956l l.l The contribution 
of Einstein led to the development of much of the field of stochastic processes, and to the notion 
that Brownian movement is due to the summated effect of a very large number of tiny impulsive 
forces delivered to the macroscopic particle being observed. This was also one of the most definitive 
arguments of the time for an atomistic picture of the microscopic world. 

Other ingenious experimentalists used this conceptual idea to e xplore the mac roscopic effects 



who devised 
Mazol tm'j ] for an 



of microscopic influences. One of the more interesting is due to Kappleil 
an experiment in which a small mirror was suspended by a quartz fiber (c.f 
analysis of this experimental setup). Any rotational movement of the mirror would tend to be 
counterbalanced by a restoring torsional force due to the quartz fiber. The position of the mirror 
was monitored by shining a light on it and recording the refiected image some distance away (so 
small changes in the rotational position of the mirror were magnified). Air molecules striking 
the mirror caused a transient defiection that could thus be monitored, and the frequencv of these 
collisions was controlled by changing the air pressure. Figure ITTTl taken from lKapplen (|193lh . shows 
two sets of data taken using this arrangement and offers a vivid depiction of the macroscopic effects 
of microscopic influences. 

In trying to understand theoretically the basis for complicated and irreversible experimental 
observations, a number of physicists have supplemented the reversible laws of physics with various 
hypotheses about the irregularity of the physical world. One of the first o f these, and arguab ly one 
of the most well known, is the so-called "molecular chaos" hypothesis of iBoltzmannI (jlQQ.'j V This 
hypothesis, which postulated a lack of correlation between the movement of molecules in a small 
collision volume, allowed the derivation of the Boltzmann equation from the Liouville equation and 
led to the celebrated H theorem. The origin of th e loss o f correlation s was never specified. In an 
effort to understand the nature of turbulent flow, R.uellel (|l978l . 19791 . Il98d ) postulated a type of 
mixing dynamics to be necessary. More recently, several authors have made chaotic hypotheses 
about the nature of dynamics at the mi croscopic l evel. The most prominent of these is Gallavottil 
(Il999l l. and virtually the entire book of IPorfmanl (jl999l ) is predicated on the implicit assumption 
that microscopic dynamics have a chaotic (loosely defined, but usually taken to be mixing) nature. 
All of these hypotheses have been made in spite of the fact that none of the microscopic dynamics 
that we write down in physics actually display such properties. 

Others have taken this suggestion (chaotic hypothesis) quite seriously, and attempted an ex - 
perimental confirmation. Figure [T31 shows a portion of the data, taken from Gaspard et al. ( 1998h . 
that was obtained in an examination of a microscopic system for the presence of chaotic behavior. 
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Figure 1.1: In the upper panel is shown a recording of the movement of the mirror in the iKappler 
( 193ll ) experiment over a period of about 30 minutes at atmospheric pressure (760 mm Hg). The 
bottom panel s hows the same experiment at a pressure of 4 x 10"'^ mm Hg. Both figures are from 
Kappler (|l93lh . See the text for more detail. 
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Figure 1.2: The data shown here, taken from Gaspard et al. ( 1998h . show the position of a 2.5 /im 
Brownian part i cle in water over a 300 s econd period with a samphng interval of ^ sec (see 



Gaspard et al.1 ( 1998h 



Briggs et al.l (|200lh for the experimental details). The inset figure shows 
the power spectrum, which displays a typical decay (for Brownian motion) with a;~^. 



Their data analysis showed a positive lower bound on the sum of Lyapunov exponents of the system 
composed of a macroscopic Brownian particle and the surrounding fluid. From their analysis, they 
argued t hat the Brown i an rn otion was due to (or the signature of) determinis tic microscop ic chaos. 
However, Briggs et al. ( 200lh were more cautious in their interpretation, and Mazol ( 20021 . Chapter 
18) has explored the possible interpretations of experiments like these in some detail. 

If true, the existence of deterministic chaos (whatever that means) would be an intriguing 
possibility since, if generated by a non-invertible dynamics (semi-dynamical system), it could serve 
as an explanation of a host of unresolved problems in the sciences. Most notably, it could serve 
as an explanation for the manifest irreversibility of our physical and biological world in the face of 
physical laws that fail to encompass irreversibility without the most incredulous of assumptions. In 
particular, it would clarify the foun dations of irreversible statistical mechanics, e.g . the operation of 
the se cond law of thermodynamics (lDorfmanLll999l : lGanavotI3 . ll999l : lM^i!A^ . ll989l . ll992l ^ 
199?! ). and the implications of the second law for the physical and biological sciences. 

In this paper, we have a rather more modest goal. We address a different facet of this chaotic 
hypothesis by studying how and when the characteristics of Brownian motion can be reproduced 
by deterministic systems. To motivate this, in Figures 11.31 through 11.61 we show the position (x) 
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and velocity (v) of a particle of mass m whose dynamics are described by 



dv 

m— = -Tv + J^it). (1.2) 

In Equations 11.11 and 11.21 JP" is a fluctuating "force" consisting of a sequence of delta-function like 
impulses given by 

oo 

J^{t) = mK^(,{t)5{t-nT), (1.3) 

n=0 

and ^ is a "highly chaotic" (exact, see Section deterministic variable generated by (,t+T = ^(Ct) 
where T is the hat map on [—1, 1] defined by: 

™"l2(i-y) for ye [0,1). ('-^^ 

In this paper we examine the behavior of systems described by equations like (|1.1|) through (|1.4|) 
and establish, analytically, the eventual limiting behavior of ensembles. In particular, we address 
the question of how Brownian-like motion can arise from a purely deterministic dynamics. We do 
this by studying the dynamics from a statistical, or ergodic theory, standpoint. 

The outline of the paper is as follows. Section [2l gives required background and new material 
including the definitions of a heirarchy of chaotic behaviours (ergodic, mixing, and exact) with a 
discussion of their different behaviours in terms of the evolution of densities under the action of 
transfer operators such as the Probenius-Perron operators. We then go on to treat Central Limit 
theorems and Functional Central Limit Theorems for non-invertible dynamical systems. Section 
returns to the specific problem that Equations 11.11 through 11.41 illustrate. We show how the par- 
ticle velocity distribution may converge and how the particle position may become asymptotically 
Gaussian, but for a more general class of maps than given by (|1.4j) . In Section 0] we illustrate the 
application of the results from Section |21 using a specific chaotic map (the dyadic map. Equation 
II. 4|) to act as a surrogate noise source. Section El considers the question when one can obtain Gaus- 
sian processes by studying appropriate scaling limits of the velocity and position variables, and 
the convergence of the velocity process to an Ornstein-Uhlenbeck process as the interval r between 
chaotic perturbations approaches 0. The paper concludes with a brief discussion in Section |21 The 
Appendix collects and extends general central limit theorems from probability theory that are used 
in the main results of the paper. 

2 Semi-dynamical systems 

We are going to examine the behavior illustrated in Section fusing techniques from ergodic theory, 
and closely related concepts from probability theory, applied to the dynamics of semi-dynamical 
(non-invertible) systems. In this section we colle ct together the necessarv machinery to do so. 
Much of this background material can be found in Lasota and Mackev ( 1994[ ). 



2.1 Density evolution operators 

Let {Yi,Bi, I'l) and {Y2,B2, 1^2) be two u-finite measure spaces and let the transformation T : Yi ^ 
Y2 be measurable, i.e. T-^{B2) C Bi where T^^{B2) = {T~^{A) : A G ^2}- Then we say that 
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Figure 1.3: The top panel shows the simulated position of a particle obeying Eauations ll.il through 
11.41 using Equation 13.281 while the bottom panel shows the velocity of the same particle computed 
with Eauation l3.211 The parameters used were: 7 = T /m = 10, k = 1, and r = — ^ ln(9 x 10~^) ~ 
0.932 so A = e~^'^ = 9 x 10~^. The initial condition on the hat map given by Equation 11.41 was 
yo = 0.12562568. 
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Figure 1.6: As in Figure [L3l with the parameters of Figure 1131 and an initial condition on the hat 
map (dH) of yo = 0.1678549321. 
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T is nonsingular (with respect to z/i and U2) if ^liT "^{A)) = for all A £ B2 with i^2iA) = 0. 
Associated with the transformation T we have the Koopman operator Ut defined by 

Urg = goT 

for every measurable function g : Y2 ^ R. We define the transfer operator Px : L^{Yi,Bi,i'i) 

{Y2 ,82,1^2) as follows. For any f G {Yi , ^Si , z/i ) , there is a unique element P^/ in {Y2 ,82,1^2) 
such that 

/ PTf{y)v2{dy) = [ f{y)iyi{dy). (2.1) 

Equation 12.11 simply gives an implicit relation between an initial density of states (/) and that 
density after the action of the map T, i.e. Prf- The Koopman operator Ut ■ L°°{Y2, 82,1^2) 
L°°{Yi,Bi,i'i) and the transfer operator Pt are adjoint, so 



g{y)PTf{y)u2{dy) = / f{y)UTg{y)Mdy) 

for g e L^{Y2, 82,1^2), f L\Y^,Bi,ui). 

In some special cases Equation 12.11 allows us to obtain an explicit form for Pt- Let Y2 = R, 
B2 = B^R) be the Borel fi-algebra, and 1^2 be the Lebegue measure. Let Yi be an interval [a,b] 
on the real line R, Bi = [a, b] n fi(R) and ui be the Lebesgue measure restricted to [o, b]. We will 
simply write L^{[a, b]) when the underlying measure is the Lebesgue measure. 

The transformation T : [a, 6] — > R is called piecewise monotonic if 

(i) there is a partition a = < ai < ... < ai = 6 of [a, 6] such that for each integer i = 1, 
the restriction of T to (aj_i, Oj) has a extension to [flj-i, ai] and 

(ii) |T'(2;)| > for X G (flj-i, aj, i = 1, 

If a transformation T : [a, 6] — > R is piecewise monotonic, then for / G L^{[a, b]) we have 

PTi[y) = 2^ |T'(r-i(y))| "'''"'''''^^ 

where T^^ is the inverse function for the restriction of T to (aj_i, Oj). Note that we have equivalently 

xeT-i({j/}) ' ^ ^' 

Of course these formulas hold almost everywhere with respect to the Lebesque measure. 

Let (Y, B) be a measurable space and let T : y — > y be a measurable transformation. The 
definition of the transfer operator for T depends on a given cr-finite measure on B, which in turn 
gives rise to different operators for different underlying measures on B. If is a probability measure 
on B which is invariant for T, i.e. i'{T~^(A)) = v{A) for all A £ B, then T is nonsingular. The 
transfer operator Pt : L^{Y,B,i') — > L^{Y,B,i') is well defined and when we want to emphasize 
that the underlying measure u in the transfer operator is invariant under the transformation T 
we will write T't,u- The Koopman operator Ut is also well defined for / G L^{Y.,B,u) and is an 
isometry of L^{Y,B,v) into L^{Y,B,v), i.e. ||C/t/||i = ||/||i for all / G L^{Y,B,u). The following 
relation holds between the operators Ut,Pt,u '■ L^(X-,8,u) L^{Y,B,i') 

VT,uUTf = f and UTVT,uf = E{f\T"\B)) (2.3) 
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for / G L^{Y,B,u), where E{-\T^^{B)) : L^{Y,B,u) L\Y,T~^{B),i^) denotes the operator of 
conditional e xpectation (see A ppendix). Both of these equations are based on the following change 
of variables dBiningslevl . Il99,4 Theorem 16.13): / G L^(y, B, v) if and only if / o T G L^{Y, B, u), 



in which case the following holds 

/ foTiyMdy)= [ fivHdy), AgB. (2.4) 
Jt-i(a) Ja 

If the measure j> is finite, we have L^iY^B^u) C L^{Y,B,i') for p > 1. The operator Ut '■ 
LP{Y,B,u) L^{Y,B,i') is also an isometry in this case. Note that if the conditional expec- 
tation operator E{-\T-'^{B)) : L^{Y,B,v) L\Y,B,u) is restricted to L'^{Y,B,iy), then this is the 
orthogonal projection of L'^{Y,B,u) onto L'^{Y,T^^{B),u). 

One can also consider any cr-finite measure m on B with respect to which T is nonsingular and 
the corresponding transfer operator Pt : L^{Y,B,m) — > L^{Y,B,m). To be specific, let 1" be a 
Borel subset of with Lebesque measure m and B = B{Y) be the cj— algebra of Borel subsets of 
Y. Throughout this paper m will denote Lebesque measure and (Y) will denote {Y, B, m). The 
transfer operator Pt : L^(Y) — > L^iX) is usually known as the Probenius- Perron operator. A 
measure v (on Y) is said to have a density g^, if i'{A) = g^:{y)dy for all A ^ B, where G L^{Y) 
is nonnegative and Jy g*{y)dy = 1. A measure v is called absolutely continuous if it has a density. If 
the Frobenius- Perron operator Px has a nontrivial fixed point in L^{Y), i.e. the equation Pt/ = / 
has a nonzero solution in L^{Y), then the transformation T has an absolutely continuous invariant 
measure u, its density 5* is a fixed point of Pt, and we call g^, an invariant density under the 
transformation T. The following relation holds between the operators Pt and 'Pt,u 

Prifg*) = 9*VT,uf for / G L\Y,B,u). (2.5) 

In particular, if the density 5* is strictly positive, i.e. g*{y) > a.e. for y £ Y, then the measures 
m and v are equivalent and we also have 



PtU) = 9*Vt,u for / G L\Y). 



The notion of a piecewise monotonic transformation on an interval can be extended to "piecewise 
smooth" transformations T : Y ^ Y with Y C R^. Therefore if T has, for example, a finitely 
many inverse branches and the Jacobian matrix DT[x) of T at x exists and det DT{x) 7^ for 
almost every x, then the Probenius-Perron operator is given by 

for / G L^iY). If T is invertible then we have Pt/(?/) = f{T~^{y))\ det DT~^{y)\. 

Finally, we briefiy mention Ruelle's transfer operator. Let y be a compact metric space, T : 
y — > y be a continuous map such that T~^({y}) is finite for each y £ Y and let : y ^ R be 
a function (typically continuous or Holder continuous). The so called Ruelle operator acts on 
functions rather than on L^{Y) elements and is defined by 

-c^/Cy) = E e'^^-V(x) 

for every y G y. The function is a so called potential. If we take ijj{x) = — log | det Dr(x)|, 
provided it makes sense, then we arrive at the representation for the Frobenius-Perron operator. 
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2.2 Probabilistic and ergodic properties of density evolution 

Let {Y, B, u) be a normalized measure space and let T : 1" — > 1" be a measurable transformation 
preserving the measure v. We can discuss the ergodic properties of T in terms of the convergence 
behavior of its transfer operator Vt.u '■ L^O^jB,!^) L^{Y,B,i'). To this end, we note that the 
transformation T is 

(i) Ergodic (with respect to z^) if and only if every invariant set A £ B is such that = or 

\ A) = 0. This is equivalent to: T is ergodic (with respect to i') if and only if for each 
/ G L^{y, B, v) the sequence ^ l^fc=o uf weakly convergent in L^{Y, B, v) to / f{y)i>{dy), 
i.e. for all g G L'^{Y,B,u) 

n—l 

i™o / llY.'PT,J(y)9iyMdy) = I f{y)y{dy) J g{y)v{dy)- 

k=0 

(ii) Mixing (with respect to v) if and only if 

lim u{A n T-''\B)) = u{A)u{B) for A,B eB. 

n—foo 

Mixing is equivalent to: For each / £ L^{Y,B,i') the sequence P^f is weakly convergent in 
L\Y,B,u) to J f{y)u{dy), i.e. 

\\TnfvlJ[y)g[y)u{dy)=ff[y)u{dy)fg[y)v{dy) for g£L^[Y,B,u). 

n^ooj J J 

(iii) Exact (with respect to i') if and only if 

lim i/(r"(^)) = 1 for ^ G 5 with T{A) G B, u{A) > 0. 

n— >oo 

Exactness is equivalent to: For each / G L^{Y,B, v) the sequence Pj- i,f is strongly convergent 
m L\Y,B,v) to J f{y)v{dy), I.e. 

lim / \V^^J{y)- [ f{y)uidyMdy) = Q. 

n^oo J J 

The characterization of the ergodic properties of transformations through the properties of the 
evolution of densities requires that we know an invariant measure i' for T. Examples of ergodic, 
mixing, and exact transformations are given in the following. 

Example 1 The transformation on [0, 1] 

T{y) = y + 4> mod 1, 

known as rotation on the circle, is ergodic with respect to the Lebesgue measure when (p is irrational. 
The associated Frobenius-Perron operator is given by 

Prfiy) = f{y - <P). 
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Example 2 The baker map on [0, 1] x [0, 1] 



0<z<i 
(2y-l,i + iz) \<z<l 

is mixing with respect to the Lebesgue measure. The Frobenius-Perron operator is given by 



/(Jy,2z) o<z<i 

/(| + iy,2z-l) \<z<l 



Example 3 The hat map on [— 1 , 1] defined by Equation 11.41 is exact with respect to the Lebesgue 
measure, and has a Frobenius-Perron operator given by 



Prfiy) 



Example 4 A class of piecewise linear transformations on [0, 1] are given by 



TN{y) 



for y £ 
for y G 



2n 2n+l' 



N ' N 
N ' N 



2n+l 2n+2~j 



(2.6) 



where n = 0, 1, . . . , [(A^ — l)/2] and [z] denotes the integer part of z. For N > 2, these piecewise 
linear maps generalize the hat map, are exact with respect to the Lebesgue measure, and have the 
invariant density 

5*(y) = i[o,i](y)- (2.7) 



Example 5 The C hebv shev maps JAdler and RivlinI ( 19641 ^) on [—1, 1] studied by Beck and Roepstorf 
(jl 987l V iReckl ()l99(il ^ and iHilgers a,nd Beckl (Il999l ^ are given by 



SN{y) = cos(A arccos y). 



A = 0,1, 



(2i 



with 5*0 (y) = 1 and Si{y) = y. They are conjugate to the transformation of Example |11 and satisfy 
the recurrence relation Sj\[+i{y) = ^yS^iy) — Siy-i{y)- For N > 2 they are exact with respect to 
the measure with the density 

9*{y) = — , ^ - 

TTy^l - y^ 

For N = 2 the Frobenius-Perron operator is given by 



PsJiy) 



2V2yT2 



/ 



2^+2 



and the transfer operator by 

Ps„uf{y) 



2^+2 



The construction of a trans fer operator for a conjugate map is presented in the next theorem 
from Lasota and Mackev ( 19941 ^ . 



13 



Theorem 1 \Lasota and Mackeil \l99A . Theorem 6.5.2)) Let T : [0,1] [0,1] he a measurable 
and nonsingular (with respect to the Lebesgue measure) transformation. Let v : i3([a, 6]) — > [0, oo) 
be a probability measure with a strictly positive density (7*, that is g^:{y) > a.e. Let a second 
transformation S : [a, b] — > [a, b] be given by S = oT o G, where 

G{x) = I g^{y) dy, a < x <b. 

J a 

Then the transfer operator Vs,u is given by 

Vs,uf = UcPTUc-^f, for f G L\[a,b],B{[a,b]),iy), (2.9) 

where Ug,Uq~i are Koopman operators for G and G~^ , respectively, and Pt is the Frohenius- 
Perron operator for T. As a consequence, v is invariant for S if and only if the Lebesgue measure 
is invariant for T. 

Example 6 The dyadic map on [—1, 1] is given by 



and has the uniform invariant density 



9*{y) = 



Like the hat map, it is exact with respect to the normahzed Lebesgue measure on [—1,1]. It has a 
Frobenius-Perron operator given by 



Example 7 [Alexander and Yorkj ( 19841 ^ defined a generalized baker transformation (also known 
as a fat/skiny baker transformation) S*^ : [—1,1] x [—1,1] — > [—1, 1] x [—1,1] by 

Spix,y) = iPx + {l-P)h{y),T{y)) 

where < /? < 1, T is the dyadic map on [—1, 1], and 



Hy) 



1, y>o, 

-1, y<0. 



For every /3 G (0,1) the transformation Sp has an invariant probability measure on [—1,1] x 
[—1,1] and is mixing. The invariant measure is the product of a so called infinitely convolved 
Bernoulli measure (see Section 0]) and the normalized Lebesgue measure on [—1,1]. If (3 = ^, 
the transformation is conjugated through a linear transform of the plane to the baker map of 
Example 121 If /3 < ^, the transformation Sfs does not have an invariant density (with respect to 
the planar Lebesgue measure). 

Example 8 The continued fraction map 

T(y) = i modi yG(0,l] (2.11) 
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has an invariant density 



g*{y) 



1 



(l + y)ln2 

and is exact. The Frobenius- Perron operator is given by 

_ 1 



(2.12) 



and the transfer operator by 



VT,uf{y)=Y.- 



y + l 



1 



{y + k){y + k + iy \y + k 

Example 9 The quadratic map is given by 

Tp{y) = l-py^, ye [-1,1] 

where < /3 < 2. It is known that there exists a positive Lebesgue measure set of parameter 
values P such that th e map h as ari abso lutely continuous (with respec t to Lebesgue measure) 
invariant measure up i .TakobsonI (|l98ll ) and Benedics and CarlesonI (jlQSfJ )). Let a > be a very 
small number and let 

A, = {/? G [2 - e,2] : |T^"(0)| > e""" and |(T^")'(r^(0))| > (1.9)" Vn > 0} 

for e > 0. lYound (|l992l ) proved that for sufficiently small e and for every /3 S the transformation 



Tp is exact with respect to Uj3 and this measure is supported on [r|(0), r;3(0)]. 

Example 10 The Manneville-Pomeau map Tp : [0, 1] — > [0, 1] is given by 

Tpiv) = y + y^^f^ modi, 

where (3 G (0, 1). The map has an absolutely continuous invariant probability measure up with 
density satisfying 

^<5*(y)<^ 



for some constants C2 > ci > (cf. iThaleil ()l98nl )). and is exact. 



Finally, we discuss the notion of Sinai-Ruelle-Bowen measure or SRB measure of T which 
was first conceived in the setting of Axion i A diffeomorphisms on con i pact Riemannian manifolds . 
This notion varies from author to author ([Alexander and Yorkj . 19841 : lEckmann and Ruell^ . 19851 : 
Tsuii I Il99fil : lYoun"5 120021 : Im^t et al.l . EoHi T Let {Y, p) be a compact metric space with a reference 
measure m, e.g. a compact subset of R'^ and m the Lebesque or a compact Riemannian manifold 
and m the Riemannian measure on y. If T : y — > y is a continuous map, then by the Bogolyubov- 
Krylov theorem there always exists at least one invariant probability measure for T. When there 
is more than one measure, the question arises which invariant measure is "interesting", and has 
led to attempts to give a good definition of "physically" relevant invariant measures. Though this 
seems to be a rather vague and poorly defined concept, loosely speaking one would expect that one 
criteria for a physically relevant invariant measure would be whether or not it was observable in 
the context of some laboratory or numerical experiment. 
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An invariant measure v for T is called a natural or physical measure if there is a positive Lebesque 
measure set Yq gY such that for every y gYq and for every continuous observable / : y ^ R 

lim -J2f{r{y))= / (2.13) 

i=0 

In other words the average of / along the trajectory of T starting in Yq is equal to the average of 
/ over the space Y . We can also say that for each y € Iq the measures ^ Y^Z^ ^T^iv) weakly 
convergent to v 

^ n— 1 

(see the next section for this notation). 

Observe that if ly is ergodic then from the individual Birkhoff ergodic theorem it follows that 
for every / E L^{Y^B,u) Condition 12.131 holds for almost all y G y, i.e. except for a subset of 
y of (z/) measure zero. Thus, if T has an ergodic absolutely continuous invariant measure v with 
density g^, then every continuous function / is integrable with respect v and Condition 12.131 holds 
for almost every point from the set {y G y : g*{y) > 0}, i.e. except for a subset of Lebesque measure 
zero. Therefore such ly is a physical measure for T. Not only absolutely continuous measures are 
phvsical measures. C o nside r, for example, the generalized baker transformation Sjs of Example d 
Alexander and Yorkj ( 1984h showed that there is a unique physical measure iz/s for each /3 G (0, 1). 



This measure is mixing and hence ergodic. Although for (i > \ the transformation expands areas, 
the measure might not be absolutely continuous for certain values of the parameter (3 (e.g. 

[i = ^^~~) in which case the Birkhoff ergodic theorem only implies Condition 12.131 on a zero 

Lebesque measure set. Therefore a completely different argument was needed in the proof of the 
physical property. 

In the context of smooth invertible maps having an Axiom A attractor the existence of a 
uniqu e physical measure on the attrac t or wa s firs t proved for A nosov diffeomorhisms by ISinail 
and later generalized by iRuellel dlSZi) and iBowenI Roughly speaking, these are 



maps having uniformly expanding and contracting directions and their physical invariant measures 
have densities with respect to the Lebesque measure in the expanding directions (being usually 
singular in the contracting directions). This property lead then to the characterization of a Sinai- 
Ruelle -Bowen measure. In a recent attempt to go beyond maps having an Axiom A attractor, YoungI 



(2003) additionally requires that T has a positive Lyapunow exponent a.e. The precise definition 



strongly relies on the smoothness and invertibility of the map T. Note that the generalized baker 
transformation Sp has Lyapunov exponents equal to In 2 and In (3 and the measure vp is absolutely 
continuous along all vertical directions. 

2.3 Brownian motion from deterministic perturbations 
2.3.1 Central limit theorems 



We follow the terminology of iBillingslevI (|1968l ). If ^ is a measurable mapping from a probabil- 
ity space (0,.?^, Pr) into a measurable space {Z,A), we call C a .^-valued random variable. The 
distribution of C is the normalized measure fi = Pro^^^ on {Z,A), i.e. 

fi{A) = Pt{C\A)) = Fr{uj : ^{u) £ A} = Pt{^ € A}, Ae A. 
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If Z = R^, we also have the associated distribution function of ( or ft, defined by 



F{x) = fj,{y ■.y<x} = Pr{C <x}, x £ R'', 

where {y : y < x} = {y : yi < Xi,i = 1, . . . , k} for x = (xi, . . . , rr^). The random variables ^ and ^ 
are, by definition, (statistically) independent if 

Pr{C eA,^eB} = Pr{C e A} Pv{^ g B}, 

i.e. the distribution of the pair is the product of the distribution of Q with that of ^. 

Let (Z, p) be a metric space and B{Z) be the a-algebra of Borel subsets of Z. A sequence {pn) 
of normalized measures on {Z,B{Z)) is said to converge weakly to a normalized measure if 

lim / f{z)iin{dz) = / f{z)n{dz) 



n—>oo 



Z 



for every continuous bounded function / : Z — > R. Note that the integrals f{z)n{dz) completely 
determine thus the sequence cannot converge weakly to two different limits. Note also that 
weak convergence depends only on the topology of Z, not on the specific metric that generates it; 
thus two equivalent metrics give rise to the same notion of weak convergence. If we have a family 
{fir : r > 0} of normalized measures instead of a sequence, we can also speak of weak convergence 
of Hr to /X when r goes to oo or some finite value tq in a continuous manner. This then means that 
Ht converges weakly to /i as r — > tq if and only if Pm converges weakly to /i for each sequence (r^) 
such that T„ — > Tq as n — > oo. 

If Z = R'^ and F and F„ are, respectively, the distribution functions of fi and fin, then (^fin) 
converges weakly to fi if and only if 

lim Fn{z) = F{z) at continuity points z of F. 

n—foc 

The characteristic function ip^^ of a normalized measure fi on R*"' is defined by 

ft^ii^) = / exp(i <r,z >)fi{dz), 



where i = V— 1 and < r, z >= ^^j=ifjZj denotes the inner product in R'^. The continuity 



theorem (|Billingslevl . Il96a . Theorem 7.6) gives us the following: converges weakly to fi if and 
only if 

lim (Pfi„{r) = ^fi{r) for each r G R*^. 



A sequence Cn of Z-valued random variables converges in distribution, or weakly, to a normalized 
measure fi on {Z,B{Z)), if the corresponding distributions of converge weakly to /i. This is 
denoted by 

Cn /i. 

If fi is the distribution of a random variable C, we write Cn — ^■'^ C- Note that the underlying 
probability spaces for the random variables C) Ci? C2--- may be all distinct. 

A sequence Cn of Z-valued random variables converges in probability to a Z-valued random 
variable C if 

lim Pr(p(Cn, C) > e) = for all e > 0. (2.14) 



This is denoted by 



Cn C. 
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Here all the random variables are defined on the same probability space. Note that if Condition 
Em holds then 

lim Pro(/9(Cn, C) > e) =0 for all e > 

n— >oo 

for every probability measure Pro on which is absolutely continuous with respect to Pr. In 

other words convergence in probability is preserved b y an absolu t ely co ntinuous change of measure. 
We will also frequently use the following result from iBillingslev (|l968l . Theorem 4.1): If {Z,p) is a 
separable metric space, and 

if Cn C and p{Cn, Cn) 0, then Cn -^'^ C- (2.15) 

We will write A^(0,cr^) for either a real- valued random variable which is Gaussian distributed 
with mean and variance o"^, or the measure on (R, ;B(R)) with density 

^ expf-I^V (2.16) 



Since aN{Q, 1) = Af(0, cj^) when a > 0, we can always write cjA^(0, 1) for CJ > 0, which in the case 
cr = reduces to 0. The characteristic function of A^(0, 1) is of the form (j){r) = exp (— |r^), r G R. 

Let (Ci)i>i be a sequence of real- valued zero- mean random variables with finite variance. If 
there is a > such that 



n 



f7iV(0,l), (2.17) 



then {Cj)j>i is said to satisfy the central limit theorem (CLT). Note that if the random variables 
C,j are defined on the same probability space (17, JF, Pr), we have the equivalent formulation of 1)2. 17() 
in the case of cj > 

lim Pr I w G : ^i=^ ^^"^^ < z \ = $(z), z G R, 



oo 



where 

^iz) = - 

/27r 

is the standard Gaussian distribution function. In the case of a 



$(z) = ^= / exp (-]:t^ ) dt (2.18) 



n 

Let {w{t),t G [0,oo)} be a standard Wiener process (Brownian motion), i.e. {w{t).,t G [0, oo)} 
is a family of real-valued random variables on some probability space (fi,^, Pr) satisfying the 
following properties: 

(i) the process starts at zero: w{Q) = a.e; 

(ii) for < ti < t2 the random variable w{t2) — w{ti) is Gaussian distributed with mean and 
variance t2 — ^i; 

(iii) for times ti < t2 < . . . < tn the increments w{t2) — w{ti), . . . , w{tn) — w{tn-i) are independent 
random variables. 
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Existence of the Wiener process {'w{t),t G [0, 1]} is equivalent to the existence of the Wiener measure 
W on the space C[0, 1] of continuous functions on [0, 1] with uniform convergence, in a topology 
which makes C[0, 1] a complete separable metric space. Then, simply, W is the distribution of a 
random variable W : Q ^ C[0, 1] defined by W{lo) : 1 1— > w{t){Lo). 

Let D[0, 1] be the space of right continuous real valued functions on [0, 1] with left-hand limits. 
We endow D[0, 1] with the Skorohod topology which is defined by the metric 



ps(V',V') = mf sup \^{t) - ^{s{t))\ + sup \t-s{t)\], € D[0,1], 

ytG[o,i] te[o,i] J 

where S is the family of strictly increa sing, continuous mappings s of [0,1] onto itself such that 



s(0) = and s(l) = 1 (jBillingslevl . 1 1 968l . Section 14). The metric space {D[0, 1], ps) is separable and 
is not complete, but there is an equivalent metric on D[0, 1] which turns D[0, 1] with the Skorohod 
topology into a complete separable metric space. Since the Skorohod topology and the uniform 
topology on C[0, 1] coincide, W can be considered as a measure on D[0, 1]. 

A stronger result than the CLT is a weak invariance principle, also called a functional central 
limit theorem (FCLT). Let {Cj)j>i be a sequence of real- valued zero-mean random variables 
with finite variance. Let a > and define the process {ipn{t),t G [0, 1]} by 

^ [nt] 

V'„(t) = ^J^O for tG [0,1], 

(where the sum from 1 to is set equal to 0). Note that ^„ is a right continuous step function, a 
random variable of D[0, 1] and ipn{0) = 0. If 

(here the convergence in distribution is in D[0, 1]), then (Cj)j>i is said to satisfy the FCLT. 

If for every k > 1 and every vector (ti, . . . , t^) G [0, l]'^ with ti < . . . < the joint distribution 
of the vector (ipniti), ■ ■ ■ i V'n(ifc)) converges to the joint distribution of {w{ti), . . . ,w{tk)), then we 
say that the finite dimensional distributions of ipn converge to those of w. For one dimensional 
distribution this convergence is equivalent to the central limit theorem. 

The convergence of all finite-dimensional distributions of ipn to those of w is not su f ficient to 



conclude that V'n w in D[Q, 1]. According to Theorems 15.1 and 15.5 of iBillingslevI (|l968l ) if, 
additionally, for each positive e 

limlimsupPr( sup \'il)n{s) - iin{t)\ > e) = {), (2.19) 

5^0 n^oo \t-s\<5 

then ^„ converges in distribution to the Wiener process w. 

The term functional central limit theorem comes from the mapping theorem ( Billineslev , 19681 . 



Theorem 5.1), according to which for any functional / : D[0, 1] — > R, measurable and continuous on 
a set of Wiener measure 1, the distribution of /(V'n) converges weakly to the distribution of f{w). 
This applies in particular to the functional /(V') = supo<5<i ip{s). Instead of real- valued functionals 
one can also consider mappings with values in a metric space. For example, this theorem applies 
for any / : Z)[0, 1] ^ Z)[0, 1] of the form /(,^)(i) = sup,<i (/>(s) or /(</.)(t) = (t>{s)ds. 
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2.3.2 FCLT for noninvertible maps 

How can we obtain, and in what sense, Brownian-like motion from a (semi) dynamical system? This 
question is intimately connected with Central Limit Theorems (CLT) for non-invertible systems 
and various invar iance principles. 

Many CLT results and invariance principles for maps have been proved, see e.g. the survey 
These esults extend back over som e decades, inc l uding contributions bv Ratnei] 



Denker (Il989h. These results extend back over som e decades, inc l uding contributions by Ratnei 
dlQT.-jVlBoyarsky and Scarowskv (|l97fl) JWonj (llQTfll 'l. Edi^ (|l98nl ^. ljIblohski and MalczakI (llQSS^ 'l. 
Jabloriskil (|l99lh . iLiveranil and I Vi anal (llQATf T 



First, however, remember that if we have a time series y{j) and a bounded integrable function 
h : X ^ R, then the correlation of h is defined as 

Rh{n) = lim — V h{y{j + n))h{y{j)). 

N^oo iV ^ — ' 
j=0 

If the time series is generated by a measurable transformation T :¥ ^ Y operating on a normalized 
measure space {Y, B, v), and if further u is invariant under T and T is ergodic, then we can rewrite 
the correlation as 

Rh{n) = lim - V h{y{j + n))h{y{j)) = / h{y)h{T^{y))v{dy). 
The average < /i > is just 

I 

<h> 

j=0 



1 ^"^ f 



Let {Y^B^u) be a normalized measure space, and T : y — > y be a measurable transformation 
such that V is T-invariant. (Y, B, v) will serve as our probability space (fi, J^, Pr). Let h £ L'^{Y, B, v) 
be such that / h{y)v{dy) = 0. The random variables Q = ho T^~^, j > 1 are real-valued, have 
zero-mean and finite variance equal to \ \h\\2 = J h?{y)i'{dy). Thus the terminology from Section 
12 .3. II applies. The explicit formulae for the Frobenius-Perron operators in Section [2. 21 show that the 
equations P^/i = or Vx^uh = can be easily solved. For instance, in Example |31 every function h 
which is odd is a solution of these equations. In particular, considering h with VT,uh = turns out to 
be very fruitful. Statistical properties of the sequence {h o T^)j>Q are summarized in the following 

Theorem 2 (CLT) Let {Y,B,v) he a normalized measure space and T : Y ^ Y be ergodic with 
respect to v. If h £ L'^{Y,B,i') is such that VT.yh = 0, then 

(i) J h{y)u{dy) = and J h{y)h{T''{y))iy{dy) = for all n > 1. 



(ii) >'^cjiV(0,l) and a 



n 



2- 



(Hi) If a > then (ho T-')j>o satisfies the CLT and FCLT. 

(iv) If h £ L°°(y, ;B, z^) and a > then all moments of — - — ■= converge to the corresponding 



moments of aN{0, 1), i.e. for each k > 1 we have 



n— >oo 



Vra ivy/ 
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Proof. First note that the transfer operator 'Pt,u preserves the integral, i-e J VT,uh{y)y{dy) = 
J h{y)v{dy). Hence J h{y)v{dy) = 0. Now let n > 1. Since z/ is a finite measure, the Koopman and 
transfer operators are adjoint on the space L^{Y,B,u). This implies 



/ 



h{y)h{T^{y)Mdy) = / hiyp^yHdy) 



and completes the proof of (i). Part (ii) follows from Lemma ^ and Theorem 1121 in Appendix, since 
for each n > 1 we have 

^ n— 1 ^ n 

— V /i o T^' = — V /i o T"-^'. 



n — ' vn 



If cj > 0, then the CLT is a consequence of part (ii), while the FCLT follows from Lemma |21 and |31 
in Appendix. Final ly, the existence and convergence of moments follow from Theorem 5.3 and 5.4 



of iBillingslevI ()l968l ) and from Lemma ^ in Appendix. 



Remark 1 Note that if T is invertible then the equation VT,u{h) = has only a zero solution, so 
the theorem does not apply. 

We now address the question of solvability of the equation VT,v{h) = 0. 

Proposition 1 Let (y, B., u) he a normalized measure space and T : Y ^ Y be a measurable map 
preserving the measure u. Let y = Yi U ^2 with Yi,Y2 G B and i'{Yi D Y2) = and let a bijective 
map If : Yi ^ Y2 be such that both ip and ip^^ are measurable and preserve the measure u. Assume 
that for every A £ B there is B £ B such that B QYi and 

T^\A) = BUip{B). (2.20) 

U h{y) + h{ip{y)) = for almost every y G Yi, then VT,uh = 0. 

Proof. From Condition 12.201 it follows that 

/ hiyHdy) = [ h[y)y{dy) + [ h{y)u{dy). 
Jt-^{A) Jb J<p{B) 



Since 



h{y)u{dy)= [ h{ip-\^{y)))u{dy), 
JviB) 



the last integral is equal to 



/ Kv>{y)){^ ° V>){dy) = / h{ip{y))u{dy) 
JB JB 



by the change of variables applied to <f ^ and finally by the invariance of v for ip. This, with the 
definition of Pt,!/, completes the proof. 
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Remark 2 The above proposition can be easily generalized. For example, we can have 

T-\A) = BUipiiB)U^2iB) 

with the sets B , ipi{B) , ip2iB) pairwise disjoint. 

Note that ifY is an interval then it is enough to check Condition \2. 2(A for intervals of the form 
[a,h). 

Example 11 For an even transformation T on [—1, 1] with an even invariant density we can take 
Yi = [—1,0] and 99(2/) = —y. In this case VT,uh = for every odd function on [—1, 1]. In particular, 
this apphes to the tent map and to the Chebyshev maps Sn of Example El with even. We also 
have Vsn,!^^ — when h{y) = y and Sjy is the Chebyshev map with N odd. Indeed, first observe 
that by Theorem^we have 'Psf^,uh = if PTj^f = where T/v is given bv 12.61 and f{y) = cos(7ry) 
for y e [0, 1] and then note that -Ptjv/ = follows because the expression 

(7V-l)/2 ^ ^ 



n=l ^ ^ ^ 



N 



reduces to 



cos(7ry)(l + 2 ^ 

n=l ^ ^ 



which is equal to 0. For the dyadic map, we can take = y + 1. Then any function satisfying 
h{y) + h{y + 1) = gives a solution to VT,uh = 0. 

The next example shows that the assumption of ergodicity in Theorem |21 is in a sense essential. 

Example 12 Let T : [0, 1] [0, 1] be defined by 

22/, yG[0,i) 
Tiy) = { 2y-l yG[i,f 



2y-l, 



The Frobenius-Perron operator is given by 



Prfiy) = If [Iv) i[o,i)(y) + ^/ {ly + \)+lf {^y + 1) %niy)- 

Observe that the Lebesgue measure on ([0, 1], ;B([0, 1])) is invariant for T and that T is not ergodic 
since T-^{[0, ^]) = [0, i] and T-'^{[^, 1]) = 1]. Consider the following functions 



hi{y) 



{ 1, ye[o,i) 
-1, yG[|,f), h2{y) 

1, ye [1,1], 



f 1, yG[o,i) 

-1, ye[|,|), 

-2, ye[|,|), 

2, ye[|,i]. 



We see at once that Py/ii = Pt^2 = 0, Pt^i = /if, and Prh^ = h^. It is immediate that for 
every y G [0, 1] we have h\{T[y)) = h\{y) = 1 and h'^iT{y)) = h\[y). Lemma^and Theorem [T^ in 
Appendix show that 



^ 71—1 

Ar(o,i), 



i=o 
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while 

n-l 



1 ^/.^oT^-^'^C, 



3=0 

where ^ has the characteristic function of the form 



1 / ^2 

^ '2, 



(f)^{r) = I exp ( -—h^iy) ) dy. 







The density of C is equal to 



11 / 1 1 / 

Consequently, the sequence (/12 o T^)j>o does not satisfy the CLT. 

We now discuss the problem of changing the underlying probability space for the sequence 
(/i o T-^)j>o. The random variables h o in Theorem |21 are defined on the probability space 
{Y,B,u). Since v is invariant for T, they have the same distribution and constitute a stationary 
sequence. We shall show that the result (iii) of Theorem [21 remains true if the transformation T is 
exact and hoT^ are random variables on {Y^B^uq) where uq is an arbitrary normalized measure 
absolutely continuous with respect to u. In other words we can consider random variables h[T^ [^q)) 
with ^0 distributed according to uq. Now these random variables are not identically distributed 
and constitute a non-stationary sequence. For example, consider the hat map T{y) = 1 — 2\y\ on 
[— 1, 1] and h(y) = y. Then Theorem |21 applies if ^0 is uniformly distributed. We will show that we 
can also consider ^0 having a density with respect to the normalized Lebesgue measure on [—1, 1]. 

Theorem 3 Let (y, i3, u) he a normalized measure space and T : Y ^ Y be exact with respect to 
V. Let h € L'^{Y,B,u) he such that VT,uh = and a = ||/i||2 > 0. If ^0 is distributed according to a 
normalized measure uq on (Y,B) which is absolutely continuous with respect to u, then (/ioT'-'(^o))j>o 
satisfies both the CLT and FCLT. 

Proof. Let go be the density of the measure with respect to u. On the probability space (1", B) 
define the random variables by 

C^{y) = J-Y^h{P{y)), yeY 

To prove the CLT we shall use the continuity theorem ( Billingslevl . 19681 . Theorem 7.6), and to do 
so we must show that 

Jirn^ j ex.p{irCn{y))9o{y)i^{dy) = exp (^-y^ , r G R. 
Fix e > 0. Since T is exact, there exists m > 1 such that 

Define 

^ n—m—l 
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and observe that 



^ m—l 



for sufficiently large n. Since for every y the sequence ^-7^ Sj^o ^(^"'(y)) converges to as n 
we obtain 



00, 



lim / exp{irCn{y))go{y)j^{dy) - / exp{irCn{T"^{y)))go{y)j^{dy) 

n^oo J J 

The equality exp(ir^„ o T"^) = C/™(exp(ir^„)) implies that 

eMirUT"'iy)))goiyHdy) = I eMirL{y))v^,ugo{y)y{dy), 



since the operators Ut and Vt,v are adjoint. This gives 



exp{irCn{T"^{y)))go{y)i^{dy) - / exp(irC„(y))i/(dy) 



< / l^?:.5o(y) - MHdy) < e. 



From Theorem [21 and the continuity theorem ( Billineslev . 19681 . Theorem 7.6), it follows that 



lim / exp{irCri{y))T^{dy) = exp ( -y ) . 



Consequently, 



lim sup 



exp{irCn{y))go{y)'^{dy) - exp ( - y 



which leads to the desired conclusion, as e was arbitrary. 

Similar arguments as above, the multidimensional version of the continuity theorem, and Lemma 
121 in Appendix allow us to show that the finite dimensional distributions of 



j=0 



converge to the finite dimensional distributions of the Wiener process w. By LemmajSlin Appendix 
Condition 12.191 holds with Pr = z^. Since z^q is absolutely continuous with respect to ly, it is easily 
seen that this Condition also holds with Pr = z^O) which completes the proof. 

Does the CLT still hold when h does not satisfy the equation VT,uh = 0? The answer to this 
question is positive provided that h can be written as a sum of two functions in which one satisfies 
the assumptions of Th eorem B] while the other is irrelevant for the Central Limit Theorem to hold. 
This idea goes back to GordinI 



Theorem 4 Let (Y, B, u) be a normalized measure space, T : Y ^ Y be ergodic with respect to 
V, and h G L'^{Y,B,v) be such that J h{y)v{dy) = 0. If there exists h E L'^{Y,B,u) such that 
'PT,uh = and the sequence Sj=o — h) o is convergent in L'^(Y,B, v) to 0, then 



n— >oo n 



(2.211 
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and 

n-l 



Moreover, if the series "^JLi f h{y)h{T'-^ {y))v{dy) is convergent, then 

/CO p 
h\y)u{dy) + h{y)h{T^{y))v{dy). (2.22) 

n=l '' 



Proof. We have 

Y:;Z^hoT^ EUhoT^ ^ EUih-h)oT^ ^ ^^^^^^ 



n \/n \ n 



Since VT,vh = 0, we obtain J h{y)u{dy) = and / hiT"^ {y))hiT^ {y))u{dy) = for i / j by Condition 
(i) of Theorem Hence 

-. n — 1 ^ n— 1 

\\h\\l='-Y.\\^oT% = \\l=Y.'hoT% 
j=0 ^ j=0 

^ n—l 

and therefore Equation 12.211 holds. Since the sequence —= } (h — h) o P is convergent to 

L^(y, i3, i^), it is also convergent to in probability. Equation 12.231 Condition (ii) of Theorem [21 
applied to /i, and property (|2.15j) complete the proof of the first part. 
It remains to show that 

lim ^= / h\y)v{dy)^2y^ / h{y)h{T^{y))v{dy). 

n—*oo n J J 



m 



n=l • 



Since v is T-invariant, we have 

, „ /n-l \ ^ n-l j 

- / 5^/.(r^(y)) u{dy) = \\h\\2 + 2-Y,Yl / Hy)h{T^{y)Hdy), 

^ \j=0 / j=l 1=1 

but the sequence (X]^=i / ^{y)h{T-' {y))^{dy))n>i is convergent to YlY=i / h{y)h{T^ {y))u{dy)^ which 
completes the proof. 

Remark 3 Note that in the above proof of the CLT we only used the weaker condition that the 

^ n— 1 

sequence —= > ih — h) o is convergent to in probability. The stronger assumption that this 

sequence is convergent in L'^{Y,B,i') was used to derive Eauation \2.21\ Note also that all of the 
computations are useless when \ \h\\2 = and the most interesting situation is when h is nontrivial, 
i.e. \\h\\2 > 0. 



Strengthening the assumptions of the last theorem leads to the functional central limit theorem. 

Theorem 5 Let {Y,B, u) be a normalized measure space, T :¥ ^ Y be ergodic with respect to v, 
and h G L'^{Y,B,u) be such that J h{y)v{dy) = 0. If there exists a nontrivial h G Lp'iY^Bjv) such 
that VT,uh = and the sequence X]"=o — h)o is v-a.e. convergent to 0, then {h o T^)j>o 
satisfies the FCLT. 
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Proof. Since every sequence convergent v almost everywhere is convergent in probability, the CLT 
follows by the preceding Remark and Theorem 01 To derive the FCLT define 

H-i H-i 
^^(t) = ^ V hoT^ and V„(t) = ^ V /ioT^ tG[0,l], 



where a = ||/i||2- Then by (iii) of Theorem |21 we have 
By property 1)2. 15(1 it remains to show that 



Psitpn,1pn) 0. 

To this end observe that 



k-l 



ps{ipn,ipn) < sup iV'n(i) " ipn{t)\ < — ^ max | V^(/l - h) O T^\. 

0<t<l CrV"^ l<k<n ^ 

— J=0 

Since the sequence X^j=o — h) o is v-a.e. convergent to 0, the same holds for the sequence 
maxi</(.<„ I "^^Zoih — h) o by an elementary analysis, which completes the proof. 

Remark 4 With the settings and notation of Theorem^and respectively Theorem\^ ifT is exact, 
then the same conclusions hold for the sequence {h o T''(^o))j>o provided that is distributed 
according to a normalized measure vq on {Y^v) which is absolutely continuous with respect to v. 
Indeed, since convergence to zero in probability is preserved by an absolutely continuous change of 
measure, we can apply the above arguments again, with Theorem replaced by Theorem 

One situation when all assumptions of the two preceding theorems are met is described in the 
following 

Theorem 6 Let (Y, B, u) be a normalized measure space, T -.Y be ergodic with respect to v , 
and h £ L'^ {Y, B , v) . Suppose that the series 

oo 
n=0 

is convergent in L'^{Y, B, u). Define f = X^^^ '^t ^''^^ h = h + f — foT. 

Then h G L^{Y,B, v), VT,uh = 0, J h{y)u{dy) = 0, and the sequence (-^ X]j=o — h)o T^)n>i 
is convergent to both in L'^{Y,B,u) and v—a.e. 

In particular, \\h\\2 =0 if and only ifh = foT — f for some f £ L'^(Y,B, u). 

Proof. Since VT,u{h + f) = f, we have by Equation 12.31 

rT,uh = VtA^ + /) - VrAf °T) = f- VT,uUTf = 0. 
Thus it remains to study the behavior of the sequence ■^X]"=o(^ — h) o T^ , which with our 
notations reduces to -^(/ o — /). This sequence is obviously convergent to in L^{Y,B,u) 
because ||/ o T" — /II2 < 2||/||2. It is also z^-a.e. convergent to which follows from the Borel- 
Cantelli lemma and the fact that for every e > the series ^{f'^°T^ > ne) = Y^^=i ^(/^ — '^^) 

is convergent as / G L'^{Y,B,v), which completes the proof. 

Summarizing our considerations for general h we arrive at the following sufficient conditions for 
the CLT and FCLT to hold. 
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Corollary 1 Let (Y, B, u) he a normalized measure space, T : Y ^ Y be ergodic with respect to v, 
and h G L'^{Y,B,u). If 

oo 

Y,\\T't,M2<^, (2.24) 

n=0 

then a > given by 

2 



h\y)v{dy) + 2Y, / h{y)h{T^{y))v{dy) 

n=l '' 

is finite and {h o T^)j>Q satisfies the CLT and FCLT provided that a > 0. 

Proof. Since the operators Vt,u and Ut are adjoint on the space L'^{Y,B, u), we have 

h{y)h{T^{y)Mdy) = [ V^MyMyMdy). 



Thus 

' ' ' < / \V^My)h{y)Hdy) < m.hMhlU 



h{y)h{T^{y))u{dy) 
by Schwartz's inequahty. Hence 

oo „ oo 

/ \h{y)h{T-{y))\dv <\\h\\2Y,\\'Pl.h\\2, 

n=l-' n=l 

which shows that the series Yl^=i I ^{y)h{T^{y))T-'{dy) is convergent. Since assumption 12.241 im- 
plies that the series Y1'^=o'Pt absolutely convergent in L'^(Y,B,i'), the assertions follow from 
Theorems HJ El and El 

Remark 5 Note that if Condition \2.24\ holds then 

hm ||P->||2 =0. 

n— >oo ' 

Since v is finite, we have 

lim =0. 

n— too ' 

Therefore the validity of Condition \2. 24\ on a dense subset of {h G L^{Y^B,v) : J h(y)u{dy) = 0} 
implies that T is exact. 

Assume that Y is an interval [a, h] in R for some o, h. Recall that a function /i : [a, 6] — > R is 
said to be of hounded variation if 



\/ h = sup^ - h{yi)\ < oo, 



1=1 



where the supremum is taken over all finite partitions, a = yo < yi < ... < yn = b, n > 1, oi Y . 

Let y([0, 1]) denote the space of all integrable functions with bounded variation over [0, 1] such 
that h{y)dy = 0. We have 

1 

\h{y)\ <\/h for he V{[0, 1]), y G [0, 1]. (2.25) 
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Example 13 For the continued fraction map there exist s a positive con stant c < 1 such that for 
every function h of bounded variation over [0, 1] we have ( losifescul . 1992I . Corollary, p. 904) 



V V^,^h <c^\/h for all n > 1, 



where i' is Gauss's measure with density as in Example |S1 From this and Condition 12.251 it 
follows that for every h G y([0, 1]) 

1 

m,M2< sup |P^Xy)| <c"\/^- 

2^e[o,i] ' 

Consequently, Condition 12.241 is satisfied and Corollary ^ applies. 

By definition, the Frobenius-Perron operator is a linear operator from L^([0, 1]) to -L^([0,1]), 
but for sufficiently smooth piecewise monotonic maps it can be defined as a pointwise map of 
y([0, 1]) into V{[0, 1]). Since functions of bounded variation have only countably many points of 
discontinuity, redefining Pt at those points does not change its properties. If, moreover, one is 
able to give an estimate for the iterates of Pt in the bounded variation norm 



BV 



/.I 

\/h+ \hiy)\dy, 

n -'0 



then obviously one is able to estimate the norm of P^f in all L^([0, 1]) spaces. In many cases there 
exist ci, C2 > and r G (0, 1) such that 

1 

mh\\BV <cir^{\/h + C2\\h\\i), heV{[0,l]). (2.26) 



We now describe two classes of chaotic maps for which one can easily show that Condition 12.241 
holds for every h £ V{[0,1]). Consider a transformation T : [0,1] [0,1] having the following 
properties 

(i) there is a partition = < ai < ... < = 1 of [0, 1] such that for each integer i = 1, ...,1 
the restriction of T to [fli-i, a^) is continuous and convex, 

(ii) r(oi_i) = and T'{ai^i) > 0, 

(iii) T'(0) > 1. 

For such transformation the Frobenius-Perron operat or has a unique fixed poin t g^:, where g^, is 
of bounded variation and a decreasing function of y (|Lasota and Mackev . 1994). Moreover it is 
bounded from below when, for example, T([0,ai]) = [0,1]. It is known ( Jabloiiski and Malczakl . 
1983bh that the estimate in Equation 12.261 holds for the Frobenius-Perron operator Pt for trans- 



formations with these three properties. Suppose that g*{y) > for a.a. y £ [0,1]. Since 
||7^T,.||oo < ll/lloo for ah / G L°°([0,l],5,i/), we have for all /i G / G L°°([0, 1], 5, i/) 

\mM2 < \\h\mV^,My' (2.27) 
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If h is of bounded variation with J h{y)g^{y)dy = 0, then hg^ G y([0, 1]) and = 
\\P^{hg,)\\i. Thus 

\\VlM^ = 0{r^/'') 

by Equation 12.261 and Corollary nappHes. 

Let a transformation T : [0, 1] — > [0, 1] be piecewise monotonic, the function jjT^^ be of bounded 
variation over [0, 1] and infj^g^,!] I^'(y)l > 1- ^'^^ such transformations the Frobe nius-Perron opera- 
tor has a fixed point g^ and g^ is of bounded variation ( Lasota and Mackev . 19941 '). Suppose that Pt 
has a unique invariant density g^. which is strictly positive. Then the transformation T is ergodic 
and g^: is bounded from below. There exists k such that T'^ is exact, the estimate in Equation 
1 2.261 is v alid for the EV obenius-Perron operator Prpk corresponding to T*^, and the following holds 
(fjablohski et al. . 19851 ^: There exists ci > and r G (0, 1) such that 

1 .1 

\P^Jiy)\<cir^(\/f+ \fiy)\dy), ye [0,1], /Gy([0,l]); 
-^0 

Hence Corollary ^ applies to T'^ and every h of bounded variation with f h(y)g^{y)dy = 0. One 
can relax the assumption that (7* is strictly positive and have instead g* > c for a.e. y G y* = {y G 
[0, 1] : > 0}. Then the above estimate is valid for y £ and / with supp/ = {y G [0, 1] : 

f{y) / 0} C y and Corollary [T] still applies to T*^. 

We now describe how to obtain the conclusions of Corollary ^ for conjugated maps. If the 
Lebesgue measure on [0, 1] is invariant with respect to T, then Theorem ^ offers the following. 

Corollary 2 Let T : [0, 1] — > [0, 1] be a transformation for which the Lebesgue measure on [0, 1] is 
invariant and for which Equation \2.2(A holds. Let g^ be a positive function and S : [a, b] — > [a, b] be 
given by S = G^^ oT o G, where 

G{x) = / g*{y)dy, a<x<b. 

J a 

7/ /i : [a, 6] — > R is a function of bounded variation with J h{y)g^{y)dy = 0, then 

00 

Y.\\viM\2<^. 

n=0 

where \\ ■ [[2 denotes the norm in L^([a, 6], i3([a, 6]), i^) and v is the measure with density g^. 
Proof. By Theorem ^ we have 

Vs,uf = UcPrUG-if, for / G L\[a,b],B{[a,b]),u). 
The operator Ug '■ -^^([0, 1]) — > L'^{[0',b],B{[a,b]),iy) is an isometry, thus 

Since G is increasing, G~^ is a function of bounded variation, as a result Uq-iH = h o G~^ is of 
bounded variation over [0, 1] , which completes the proof. 

Finally we discuss the case of quadratic maps. We follow the formulation in Vianal ( 1997 '). 
Consider the quadratic map T^, (3 G (0,2), of Example 1^1 and assume that for the critical point 
c = there are constants Ac > 1 and a > such that Ac > e^" and 
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(i) |(T^)'(T(c))| > for every n > 1; 

(ii) \Tp{c) — c| > e"""" for every n > 1; 

(iii) Tp is topologically mixing on the interval \T'^{c),Tp{c)\, i.e. for every interval / C [—1,1] 
there exists k such that T^{I) D \r'p{c),Ti3{c)]. 

lYound (Il992l ) shows that there is a set of parameters (5 close to 2 for which conditions (i), (ii), 
(iii) hold with Ac = 1.9 and a = 10~^ and that the central limit theorem holds for {h o T^) with 
h of bounded variation. Vianal ( 199?! ^ (Section 5) proves that conditions (i), (ii), (iii) imply all 
assumptions of our Corollary ^ for the transformation and functions h of bounded variation 
with J h(y)vp{dy) = 0. 

We now use the following generalization of Theorem El which allow us to have the CLT and 
FCLT for maps with polynomial decay of correlations. 

Theorem 7 \Tvmn-Kaminsk\ \200A) ) Let {Y^B, v) he a normalized measure space, T -.Y ^ Y he 
ergodic with respect to u, and h E L'^{Y,B,u) be such that f h{y)v{dy) = 0. Suppose that 

oo n—1 

Y.n-'i\\Y,VlM\2<^- (2.28) 

n=l fc=0 

Then there exists h G L'^{Y,B, u) such that Vt,u^ = and the sequence (-^ X^j=o (^ — h)o T^)n>i 
is convergent in L'^{Y,B,u) to zero and if 



k=0 



n—1 ^ 

"^V^^M2 = ©(n") with a<- 



then (-^ Yl^=o {h — h) o T^)n>i is convergent u—a.e. to 0. 

Now we give a simple result that derives CLT and FCLT from a decay of correlations. 

Corollary 3 \Tyra,n-Ka,mAnsk\ 1200^ )) Let {Y, B, v) he a normalized measure space, T : Y ^ Y he 
ergodic with respect to v, and let h £ L°^{Y,B,u) he such that J h{y)u{dy) = 0. // there are a > 1 
and c > such that 

h{y)g{T^{y))u{dy) < ^\\g\\oo (2.29) 
for all g G L^{Y,B, u) and n>l, then o" > given by 

/oo „ 
h\yMdy)+2Y, / h{yMT-{y)Mdy) 
n=l '' 

is finite and {h o T^)j>Q satisfies the CLT and FCLT provided that a > 0. 

Only recently the FCLT was established by iPollicott and for maps such as the 

Manneville-Pomeau map of Example 1 101 and for Holder continuous functions h with f h(y)i'(dy) = 
under the hypothesis that < /3 < 5. When < /3 < ^ the CLT was proved bv lYounj (|l999h . 
where it was shown that in this case condition 12.291 holds. Thus our Corollary |31 gives both the 
CLT and the FCLT for maps satisfying the following: 
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(i) T(0) = 0, r'(0) = 1, T is increasing and piecewise and onto [0, 1], 

(ii) infe<j/<i |T'(y)| > 1 for every e > 0, 

(iii) lim^^oy'"''T"(y) /O. 

as the "tower method" of lYoungl (jl99fll ) gives us the estimate in Equation 12.291 with a = — — 1 for 
all Holder continuous h and g G L°°([0, 1], u) with the constant c dependent only on h. 

2.4 Weak convergence criteria 

Let {X, I • I) be a phase space which is either R'^ or a separable Banach space, and denote by A^i 
the space of all probability measures defined on the u-algebra B{X) of Borel subsets of X. For a 
real- valued measurable bounded function /, and ;U S A^i, we introduce the scalar product notation 

(/,^) = / f{x)n{dx). 
Jx 

One way to characterize weak convergence in A^i is to use the Fortet-Mourier metric in A^i, which 
is defined by 

c?fm(/Wi,/^2) = sup{| - {f,fJ-2) I : / G ^fm} for Hi,H2 G Mi, 

where 



TFM = {f-X^-R: sup < 1, I/Il < 1} 



and = sup^^y ^"^^n^rrp^- This defines a complete metric on A^i, and we have fin ^ fJ' weakly 



if and only if dpuilJi-n, fj) — > (cf. IPudlev (|l989h [Chapter 3]) 



We further introduce a distance on A^i by 

d{fJ-i,H2) = sup{| - (/,Ai2) I : < 1} for W,/^2 G Mi. 

This quantity is always defined, but for some measures it may be infinite. It is easy to check that 
the function d is finite for elements of the set 

M\ = {// G Ml : / \x\ij,{dx) < oo}, 
Jx 

and defines a metric on this set. Moreover, Al} is a dense subset of {Mi, dpAi) and 

dFM{.l^l,IJ'2) < d{ni,fi2)- 

Let {Y,B, v) be a normalized measure space and let i?„ : X x y ^ X be a measurable trans- 
formation for each n G N. We associate with each transformation an operator P„ : Mi Mi 
defined by 

P„/i(A)= / / lA{Rn{x,y))y{dy)fx{dx) (2.30) 
Jx Jy 

for /i G A^i, where 

, . ( 1 xeA 
^^(^) = \ x^A 
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is the indicator function of a set A. Write 

Unfix) = l^f{Rn{x,y))u{dy) 

for measurable functions / : X 
satisfy the identity < Unf,fJ- >= 
X defined by 



then Unfix) =< f,Pn5x >■ 

Remark 6 Note that if Rnix,y) does not depend on y, then Unf = U^^f where Ur^ is the Koop- 
man operator corresponding to Rn : X ^ X. The following relation holds between the Frobenius- 
Perron operator Pr„ on L^iX,BiX),m) and the operator Pn-' If ^ has a density f with respect to 
m, then PR„f is a density of PnH- 

On the other hand if Rnix,y) does not depend on x, then Unf is equal to J Uii^fiy)vidy), where 
Ur„ is the Koopman operator corresponding to Rn : Y ^ X . The operator Pn has the same value 
V o R~^ for every \i ^ M.\. 

Assume that for each n G N the transformation i?„ : X x y — > X satisfies the fohowing 
conditions: 

(Al) There exists a measurable function L„ : y ^ such that 

\Rnix,y)-Rnix,y)\<Lniy)\x-x\ for yEF. 

(A2) The series 

oo » 

^/ |i?„(0,r(y))-i?„+i(0,y)|i/((iy) 

n=l 

is convergent, where T : y ^ y is a transformation preserving the measure v. 
(A3) The integral |i?„(0, y)|i^((fy) is finite for at least one n. 

Proposition 2 Let the transformations Rn satisfy conditions (A1)-(A3). If 

lim / Lniy)vidy) = 0, 

n— >oo Jy 

then there exists a unique measure /i* G Mi such that iPnfJ-) converges weakly to for each 
measure fi £ ^Al. 

Proof. Assumptions (A2) and (A3) imply that P„(A^J) C M\. By the definition of the metric d 
we have 

diPn6o,Pn+i6o) = sup{|C/„/(0) - C/„+i/(0)| : < 1}. 
Since the transformation T preserves the measure v, we can write 

UnfiO) = fiRniO, y)vidy) = fiRniO, Tiy))vidy) 



R, for which the integral is defined. The operators Un and Pn 
--< f, PnfJ, >. Note that fi = 5x, where 6x is the point measure at 
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for any / with < 1. Hence 

\Unf{0)-Un+ifm< I \Rn{0,T{y))-Rn+i{0,y)\u{dy). 

Consequently 



Y 



d{Pndo,Pn+i5o) < / \Rn{0,Tiy)) - Rn+iiO,y)\i^{dy), 



Y 



and 



dFM{PnSo,Pn+lSo) < (i(P„(5o, -P„,+l5o) ■ 

From Condition (A2), the sequence {Pn.So) is a Cauchy sequence. Since the space {Mi, dpM) is 
complete, (Pn^o) is weakly convergent to a /i* G Mi. From (Al) it follows that 



d{PnlJLi,PnlJi2) < j Ln{y)u{dy)d{^ii, ^2) ■ 



Hence {PnfJ') is weakly convergent for each n £ Ml and has the limit ^u*. Since, for sufficiently 
large n, each operator P„ satisfies 

dFM{PnlJ'l,PnlJ'l) < dpuitJ-l, 

and the set M\ is dense in {Mi,dFM), the proof is complete. 

3 Analysis 

We now return to the original problem posed in Section^ We consider the position (x) and velocity 
(v) of a dynamical system defined by 

= v{t), (3.1) 



dt 
dv{t) 

dt 



h{v{t))+i^{t), (3.2) 



with a perturbation i] in the velocity. We assume that rj{t) consists of a series of delta-function-like 
perturbations that occur at times Iq, ti,t2,- ' These perturbations have an amplitude /i(^(t)), and 
rj{t) takes the explicit form 

00 

v{t) = K^h{mnt-tn). (3.3) 

n=0 

We assume that ^ is generated by a dynamical system that at least has an invariant measure 
for the results of Section 13.11 to hold, or is at least ergodic for the Central Limit Theorem to hold 
as in Section ISI21 

In practice, we will illustrate our results assuming that £ is the trace of a highly ch aotic semidy- 



namical system that is, indeed, even exact in the sense of iLasota and Macke\ J (II 9m (c.f. Section 



nignly cn aotic 

sJdiiai) (c.f. 



121). £ could, for example, be generated by the differential delay equation 



sf^ = -m+mt-i)), (3.4) 
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where t he nonlinearity T has the appropria te properties to generate chaotic solutions ( Mackev and Glas3 
( 1977 ): an der Heiden and Mackev (jl982h ). The parameter 5 controls the time scale for these os- 
cillations, and in the limit as 5 —> we can approximate the behavior of the solutions through a 
map of the form 

= r(e„). (3.5) 

Thus, we can think of the map T as being generated by the sampling of a chaotic continuous time 
signal ^(t) as, for example, by the taking of a Poincare section of a semi-dynamical system operating 
in a high dimensional phase space. 

Let {Y,B,u) be a normalized measure space. Let h : R'^ — > R*^, h :Y ^ R'^ be measurable 
transformations, and let {tn)n>o be an increasing sequence of real numbers. Assume that ^ : 
R+ X y — > y is such that ^{tn+i) = T{^{tn)) for n > 0, where T : 1" ^ y is a measurable 
transformation preserving the measure u. Combining (|3.2|) with (|3.3() we have 

^ = h{v{t)) + K f; hmm^ - tn) m 

n=0 

We say that v{t), t > to, is a solution of Equation 13.61 if. for each n > 0, v{t) is a solution of the 
Cauchy problem 

= b{v{t)), t G (t„,t„+i) ^^ j-^ 

v{tn) = v(t-) + Kh{C{tn)), 

where vo is an arbitrary point of R*"' and v{t^) = lim^_^^- v{t) for n > 1. 

Let vr : R+ X R*^ ^ R'' be the semigroup generated by the Cauchy problem (if 6 is a Lipschitz 
map then tt is well defined) 

< ^ = ^(^(*))' ' > (3.8) 
v{0) = Vo, 

i.e. for every vq G R'^ the unique solution of 1)3. 8(1 is given by v{t) = 7r{t,vo) for t > 0. As a result, 
the solution of (|3.6|) is given by 

v{t) = 7r{t - tn, v{tn)), for t £ [tn, Wl)> " > 0- 

After integration, for t G [tn,tn+i) we have 

rt rt j-t—tn 

:{t)-x{tn)= I v{s)ds = / iris -tn,v{tn))ds = / 'K{s,v{tn))ds. 

Consequently, the solutions of Equations 13.11 and 13.21 are given by 

rt—tn 

x{t) = x{tn)+ / Tr{s,v{tn))ds, (3.9) 
JO 

v{t) = ■K{t-tn,v{tn)), for t G t„+i), n > 0. (3.10) 

Observe that x{t) is continuous in t, while v{t) is only right continuous, with left-hand limits, and 

v{tn) = hmj^j+ v{t). 



x( 
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We are interested in the variables v{tn), Vn := v{t^), ■= S.{tn), and x„ := x{tn) which appear 
in the definition of the solution v{t) and x{t). We have 

V{tn) = Vn + Kh{^n), (3.11) 

Vn+1 = -Kitn+l -tn,v{tn)), (3.12) 

= r(C„), (3.13) 

rtn + l—tn 

Xn+l = Xn+ 'IT{s,v(tn))ds. (3.14) 



Jo 

We are going to examine the dynamics of these variables from a statistical point of view. 

Suppose that vo has a distribution has a distribution v, and that the random variables 

are independent. What can we say about the long-term behavior of the distribution of the random 
variables v(tn) or 

3.1 Weak convergence of v{tn) and Vn- 

To simplify the presentation and easily use Proposition |21 of Section 1^31 assume that the differences 
tn+i — tn do not depend on n, and that i„ = nr for n > 0. Define A : R'^ ^ by 

K{v) = 7r(T,i;), V G R'', (3.15) 

where tt describes the solutions of the unperturbed system as defined by Equation 13.81 In par- 
ticular, adding chaotic deterministic perturbations to any exponentially stable system produces a 
stochastically stable system as stated in the following 

Corollary 4 Let A : R*^ — > R*"' he a Lipschitz map with a Lipschitz constant A G (0, 1) . Let T : Y ^ 
Y be a transformation preserving the measure u, and h -.Y ^ R'^ he such that Jy \h{y)\L'{dy) < oo. 
Assume that the random variahles vq and independent and that S^q has a distribution v. 

Then v{nT) converges in distribution to a probability measure /i* on R'^ and /i* is independent of 
the distribution of the initial random variable vq. Moreover, Vn converges in distribution to the 
probability measure /i* o A^-*^. 



Proof. From Equations 13. Ill and 13. 12l it follows that 

v{{n + l)r) = A{v{nT)) + k/i(^„+i), n > 0. 
Define the transformation i2„ : R'^ x y — > R*^ recursively: 
f Ro(.v,y) = v + Kh{y), 

\ Rn+i{v,y) = AiRniv,y)) + Kh{T^+^y)), veK',yeY,n>{). ' 

Then v{nT) = i?„(fo,Co). One can easily check by induction that all assumptions of Proposition 
l^are satisfied. Thus v{nT) converges in distribution to a unique probability measure on R'^. 
Since Vn+i = A{v{nT)) and A is a continuous transformation, it follows from the definition of weak 
convergence that the distribution of Vn+i converges weakly to fi^ o A^^. 

We call the measure fi^^ the limiting measure for v{nT). Note that /i* may depend on u. 

Remark 7 Although this Corollary shows that there is a unique limiting measure, we cannot con- 
clude in general that this measure has a density absolutely continuous with respect to the Lebesgue 
measure. See Example^ and B.emark M'iA 
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3.2 The linear case in one dimension 

We now consider Eg nation 13. 21 when b{v) = —'jv and 7 > 0. In this situation, we are considering a 
frictional force hnear in the velocity, so Equations 13. II and 13.21 become 



dv{t) 



-fv{t)+Kj2hmmt-tn). (3.17) 



dt 

n=0 

To make the computations of the previous section completely transparent, multiply Equation 
13.171 bv the integrating factor exp^-jt), rearrange, and integrate from (t„ — e) to (tn+i — e); where 
< e < min„>o(tn+i - in), to give 

v(tn+i - e)e^(*"+i-^) - vitn - e)e^(*"-^) = ^Y. KH^mz - tn) dz 

n=0 -^ti-^ 

= t,e^i'--^)h{i{tn-e)) (3.18) 
Taking the lim£_»o in Eo nation 13.181 and remembering that v{t~) = Vn and ^(in) = Cn, we have 

Vn+l = XnVn + kA„/i(C„), (3.19) 

where = tn+i — and < A„ = e"'^'^" < 1. 

We simplify this formulation by taking — t„ = r > so the perturbations are assumed to 
be arriving periodically. As a consequence, A.„ = A with 

A = e"^^. (3.20) 

Then, Eq nation 13 . 1 91 b ecomes 

Vn+l = XVn + K\h{in)- (3-21) 

This result can also be arrived at from other assumptions ^. 

3.2.1 Behaviour of the velocity variable 

For a given initial vq we have, by induction, 

n-l 

= A"7;o + kA J^A"-^-^/i(C,) 

j=0 
n-l 

= A"7;o + kA^A"-^-^7i(T^(Co)). (3.24) 

j=0 



^Alternately but, as it will turn out, equivalently, we can think of the perturbations as constantly applied. In this 
case we write an Euler approximation to the derivative in Eauation l3.17l so with an integration step size of r we have 

v{t + r) = (1 - -iT)v{t) + Tnh{£,{t)). (3.22) 

Measuring time in units of r so t„+i = t„ + r we then can write this in the alternate equivalent form 

Vn+l = \Vn + ni\h{^„), (3.23) 

where, now, A = 1 — 7r and ki — ktX"^ . Again, by induction we obtain Eauation l3.24l 
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We now calculate the asymptotic behaviour of the variance of Vn when is distributed according 
to u, the invariant measure for T. Assume for simplicity that vq = 0, set cr'^ = J h'^{y)i'{dy) and 
assume that J h{y)h{T^{y))u{dy) = for n > 1. Then we have 



j vlu{dy) = k'I I^X; A"-^7i(r^ (y)) J u{dy). 
Since the sequence hoT^ is uncorrelated by our assumption 

\j=o / j=0 

Thus 

.2 



I 



viu{dy) = K'a' ( ) . (3.25) 



Since b{v) = — jv, we have Tr{t,v) = e '^^v. Equation 13.211 in conjunction with Equations 13.121 
and I3.1l| leads to Vn+i = Xv{nT) where v{nT) = Vn + Kh{S^n)- Thus 

n 

v{nT) = A"7;o + k J]] A"~J/i(rJ (^o)) (3.26) 
i=o 

Case 1: If A < 1 and is distributed according to i^, then by Corollary 0] there exists a unique 
limiting measure for v{nT) provided that h is integrable with respect to u. The sequence Vn also 
converges in distribution. Since the function A defined by Equation 13. 151 is linear, A(u) = \v, both 
sequences v{nT) and Vn are either convergent or divergent in distribution. 

What can happen if the random variable in Equation 13.261 is distributed according to a 
different measure? 

Proposition 3 Let A < 1, the transformation T he exact with respect to v, and h G L^{Y,B,u). 
If the random variable S^o is distributed according to a normalized measure vq on {Y, B) which is 
absolutely continuous with respect to v, then 

n 

Kj;A"-^7i(r-'(eo))-"/i* (3.27) 

3=0 

and /.f* does not depend on uq. 

Proof. Recall that by the continuity theorem 

n 

K5^A"-^7i(r-'(eo))-" 

0=0 

if and only if for every r G R 

lim / exp(irft; > \^~^ h{T^ {y)))vQ{dy) = / exp(irx)/i*((ix). 

"^"^ Jy f^o 
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We know that the last Equation is true when uq = u. Since for every m > 1 the sequence 



m—l 



j=0 



is convergent to as n — > oo and 



n 



j=m. j=0 



an analysis similar to that in the proof of Theorem |31 completes the demonstration. 

-vn— 1 
^j=0 



Case 2: If A = 1 we have Vn = vq + Z]?=o ^(^'' ('^o))- Since vq and are independent ran- 



dom variables, vq and nYl]^=o ^(^"'(^o)) are also independent. Hence Vn converges if and only if 
k;^"^q h(T^{S^o)) does. Moreover, if there is a limiting measure /x* for u„, then the sequence v{nT) 
converges in distribution, say to i/*, and ^.^ is a convolution of the distribution of vq and z>*. As 
a result, depends on the distribution of vq. However, if the map T and function h satisfy the 
FCLT, then 



Hence E"=o ^(^-^ (Co)) is not convergent in distribution since the density is spread on the entire 
real line. 



3.2.2 Behaviour of the position variable 

For the position variable we have, for t G [nr, (n + 1)t), 

ft /■* _ I' _ ^ 1 - e"''''^*""'^) 
x{t) - x{nT) = I v{s)ds = I e~'^^^~'^'^'v{nT)ds = f(nr). 

J TIT J TIT T 



With x(nT) = Xn we have 



1 - e-'T' , , 1 - A , , , ^ 

Xn+i — Xn = vinr) = v{nT). (3.28) 

7 7 



Summing from to n gives 



1 - A 



n 



Xn+l = xo-\ — ^ v{jT). (3.29) 

From this and Equation 13.261 we obtain 



^ .=0 



1 - A " / ^ \ 

= X0 + ^ A% + K^A^'-^/i(r(eo)) 

^ j=0 \ i=0 ) 

(l-A"+i) (1-A)k^^A^,._. 
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Changing the order of summation in the last term gives 

n j n n " 1 — \n-i+l 

j=0 j=0 i=0 j=i 4=0 



^ n , n 



i=0 i=0 

Consequently 

= xo + + - E /i(5^'(eo)) - — E A"^^/i(T^(eo)). 



In conjunction with Eouation 13.281 this gives 

(1 - A") K 



xn = x^ + ^ + - E KT\io)) - - E A'^-^/iir (eo)). (3.30) 



Next we calculate the asymptotic behavior of the variance of Xn- Assume as before that xq 
vq = 0, and that / h{y)h{T^{y))u{dy) = and cj^ = / h?{y)v{dy). We have 



2 



EMr^(eo))-E^"^'^(^^(^o) = E^(^^(^o)) + E^""^(^*(^o)) 

\j=0 i=0 / \j=0 / \i=0 / 

n n 

-2EMr*(^o))E^""^(^*(^o)) 

i=0 1=0 

Since, by assumption, the sequence /i o T* is again uncorrelated we have 

n \ 2 n n 

Yh{riy))\ u{dy) = EE/ ^(^*(2/))/i(r^(y)Mrf?/) 

/ j=0 

E / h{riy))Mdy) = {n + l)a^ 



\i=0 / 1=0 j=0 

n. 



i=0 

Analogous to the computation for the velocity variance 

/ f'EA"-^/i(r(y))'j K^^2/) = Ea2"-2. I ^(^^ 

\j=0 / 1=0 



1 _ \2n+2 



and 



n n n »i „ 

i=o i=o 

n -, _ \n+l 

EA"-^y Mr(y))2z.(dy) = ^3^cT 



i=o i=o 4=0 i=o 



2=0 

Consequently, if xq = "Wo = then 



= — ^ (^n + 1 + ^— ^ - 2^— ^ ) . (3.31) 
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Theorem 8 Let (Y, B, v) he a normalized measure space, T :¥ ^ Y be a measurable map such that 
T preserves the measure v, let a > be a constant, and h G L'^{Y, B, v) be such that J h{y)u{dy) = 0. 
Then 

EtoHT'm^d ^^^^^2^ (3.32) 



if and only if 



Proof. Assume that Condition 13.321 holds. Prom Eouation I3.3UI we obtain 

^ = ^ + <1-^^ + -^±HTm - « A"-'/.(r'(f„)). 

By assumption 

n 

Thus the result wih follow when we show that the remaining terms are convergent in probability 
to zero. The first term 

xo ^ (1 - A") vo 



n 7 -y/n 

is convergent to zero a.e. hence in probability. The sequence X^ILo '^"~*^(^*(Co)) is convergent in 

distribution and — y= ^ as n ^ oo. Consequently, the sequence 
7V^ 



n 



is convergent in probability to zero which completes the proof. The proof of the converse is analo- 
gous. 

Remark 8 Observe that if the transformation T is exact and is distributed according to a mea- 
sure absolutely continuous with respect to u, then the conclusion of Theorem\^ still holds. 

Theorem IHl generalizes the results of Chew and Tind (j2nfl2h . In Section [2.3.21 we have discussed 
when Condition 13.321 holds for a given ergodic transformation. 

Remark 9 Note that if we multiply Gaussian distributed random variable N(0, 1) by a positive 
constant c, then it becomes Gaussian distributed A^(0, c^) with variance (? . Thus if we multiply 

both sides of Eauation \3.3'^ bu —=, we obtain 



/UT r7 
So if K = yJ^ruT, as in Chew and Tind \200A ). then 
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4 Identifying the Limiting Velocity Distribution 

Let y C R be an interval, B = B{Y), and T : Y ^ Y he a transformation preserving a normalized 
measure v on (Y,B{Y)). Recall from Section [3.11 that /x^, is the limiting measure for the sequence 
of random variables (z;(nr)) starting from vq = 0, i.e. 

n 

where h : Y ^ H is a given integrable function, < A < 1, = T^{^o), and is distributed 
according to i^. 

Proposition 4 Let Y = [a, b] and /i : y — > R 6e a bounded function. Then the limiting measure 
fiif has moments of all order given by 

x^H^^{dx) = lim / v{nT)^iy{dy) (4.1) 

n— »oo J 

and the characteristic function of ^* is of the form 

\k 



i-XJ / ' \ ft P 



fc=0 

Moreover, '^^ , 7— -r]) = 1, where c = sup gy \h{y)\. 

i — A i — A 

Proof. Since h is bounded, we have 

The existence and convergence of moments now follow from Theorem 5.3 and 5.4 of BillingslevI 



1968^ . Since v{nT) has all its values in the interval — — ^] we obtain //„([— - — — ^]) 



1 where Hn is the distribution of v{nT). Convergence in distribution ( BillingslevI . tl96^ . Theor 
2.1) implies that 

limsup^„(F) < /i=i.(F) 



em 



n^oo 



for all closed sets. Therefore ^*([ -, -]) = 1. The formula for the characteristic function 

1 — A 1 — A 

is a consequence of the other statements, and the proof is complete. 

Remark 10 Note that if the characteristic function of /i* is integrable, then //* has a continuous 
and bounded density which is given by 

1 

f*{x) = TT ex.]i{—ixr)(f)^{r)dr, x G R. 
^71" J-00 

On the other hand if has a density then (j)^:{r) ^ as \r\ 00. 

Note also that if a density exists then it must be zero outside the interval ' 



1-A' 1- A' 
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Let Y = [a, b] be an interval and h{y) = y for y G y. Thus h is bounded and for this choice 
of h all moments of the corresponding limiting distribution /U* exists by Proposition^ However, 
the calculation might be quite tedious. We are going to determine the measure /i* for a specific 
example of the transformation T by using a different method. 

Let hn ■ Y ^ Hhe defined by 

n 

/in(y) = J^A"-^r(y), yeY,n>0. (4.2) 
1=0 

Then v{nT) = Khn{^o) ^^nd v„ = kA/i„,_i(^o)- Thus knowing the limiting distribution for these 
sequences is equivalent to knowing the limiting distribution for hn{(,o)- 



4.1 Dyadic map 

To give a concrete example for which much of the preceding considerations can be completely 
illustrated, consider the generalized dyadic map defined by Equation 12.101 



2y + l, ye [-1,0] 
2y-l, ye (0,1]. 



Proposition 5 Let ^,^o be random variables uniformly distributed on [—1,1]. Let (a^) be a se- 
quence of independent random variables taking values drawn from {—1,1} with equal probability. 
Assume that ^ is statistically independent of the sequence (a^). Then for every A G (0, 1) 



Proof. The random variable 

oo 

fc=i 

is uniformly distributed on [—1,1]. It is easily seen that for the dyadic map 

oo 

e. = r^(eo) = E^ for ^^1- (4-4) 

k=l 

Using this representation we obtain 

n—l n—1 n—i n—l oo 

1=0 i=0 k=l i=0 k=n-i+l 

Changing the order of summation leads to 

n / k ,n-l-k+i\ \n-l-i 0° 

El \ ^ ^ \ \ ^ A \ ^ Ok+n _ 

k=\ \t=\ / 1=0 k=l 
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Consequently 



where 



This gives 



n-l 



w„ 



A ^ 

k=l 



n 

E 

.fc=i 



at + 



A 




(4.5) 



Note that for every n we have \wn\ < 2. Therefore X^Wn is a.s. convergent to as re ^ oo. Since 
^n(Co) converges in distribution, say to Jl^,, we have hn-i{S,o)+X""Wn — and the random variables 
on the right-hand side of Equation 14.51 converge in distribution to Jl^: . Since the random variables 
ak are independent, the random variables Ylk=i X^~^(^k and ^„ are also independent for every re. 



The same is true for Y11:=q ■^'^Qifc+i ^nd ^. Moreover, + Yl'k= 
identical distributions. Thus 



y-'^au and ^ + Ylll A'^afc+i have 



1 



A 



f n-l \ 
V fc=0 / 



On the other hand X]fc=o ^'^'^fe+i ~^ Yl'k=o '^k+i almost surely as re — > oo, but this implies 
convergence in distribution. The proof is complete. 

Before stating our next result, we review some of the known properties of the random variable 
which appears in Equation 14.31 For every A G (0, 1) let 



Ca 



oo 

E 

fc=0 



A afc+i. 



(4.6) 



and let q\ be the distribution function of Caj Q\{x) = P^ICa < x} for x € R. Explicit expres- 
sions for Q\ are, in general, not known. T he measure i n duced by the distribution q\ is called an 
infinitely convolved Bernoulli measure (see Peres et al.l (|2000l ) for the historical background and 

recent advances V 

It is known (j.Tessen and Wintnei . that Q\ is continuous and it is either absolutely contin- 

uous or singular. Recall that x is a yoint of increase of p\ if p\ {x — e) < Qx{x + e) for all e > 0. 
The set of points of increase of ( Kershner and Wintneii Il935l ) is either the interval [— Yzr^' 
when A > ^ or a Cantor set K\ of zero Lebesgue measure contained in this interval when A < ^, 
Q\ is always singular f or A £ (0 ' _^) an d the Cantor set K\ satisfies K\ = {\K\ + 1) U {\K\ — 1) 
and j-ry,— ^ -^a- Wintnei ( 19351 ^ noted that qx has the uniform density p\{x) = |l[_2,2](a^) 
for A = ^ and that it is absolutely continuous for the fc th roots of ^. Thus it was suspected that 
Qx is absolutely continuous for all A G [^>1)- However, Erdod ( 1939al l showed that qx is singular 

for A = and for the reciprocal of the so called Pisot numbers in (1,2). Later lErd5sl (|l939bl ) 

showed that there is a /3 < 1 such that f or almost all A G (/?,!) the measure qx is absolutely 
continuous. Only recently, ISolomvakI (Il99,5h showed that /3 = o- 

Proposition 6 For every A G (0, 1) the density of the limiting measure fi^ of v{nT) satisfies 



f^{v) = if and only if l-yj > 



1 - A' 
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Moreover, on the interval [—j^, j^) we have 



( /'2-A„, 



1-A 

\v\ < 



<v< 



(2-A)(l-A) 



where q\ is the distribution function of Cx defined by Equation \4-(>\ 



(2-A)(l-A) 



A 



(2-A)(l-A) < < — > 



Proof. Recall that we have v{nT) = /in('^o)- By PropositionlHl the sequence (2 — A)/i„(^o) converges 
in distribution to ^ + and the random variables ^ and (^\ are statistically independent. Since ^ 
has the uniform density on [—1, 1] and q\ is continuous, the density of + Ca is given by 



- z)dpxiz) 



Since f (nr) converges in distribution to 



OO -j^ 



,+i]iz)dpx{z 

OO 

x+1 ^ ^ 

-dpx{z) = -{qx{x + 1) - px{x - 1)). 

x-1 



K. 



A 



A 

+ Ca)- Thus p^ has a density given by 



2 r-v . 2 

+ Ca)i it follows that p^ is the distribution of 



A 



Qx 



Consequently, f^{v) 



2k 

if and only if 
2- A 

^-A 



-V + 1 



QX 



V G R. 



-v + 1 



ex 



Since qx is nondecreasing, it must be constant outside the set of points on which it is increasing, 
which is contained in the interval [— jhxl- Hence gx{x) = for x < and gx{x) = 1 for 

X > jz^- Therefore if \v\ > then f^{v) = 0. If A > ^ the function gx is increasing on the 
interval [— thus f^ is positive on [— jta' T^]- Now let A < ^. Since j^, —jhx ^ ^x, we 
also have 



i l-2A 
l-A 



2A-1 1-2A 
1-A ' 1-A 



E Kx and they divide the interval [— j^] into three intervals of length 2, 
, and 2 respectively. Since the middle interval has length less than 2 and the distance between 



the points ^—^v + 1 and ^—^v — 1 is always 2, the result follows in this case as well. 



Remark 11 If X = ^ then gx{x) = i(x + 2), 2 — A 



9 
32k 
_3_ 
8k' 

9 
32k 



{v + 2), -2k<v< 



^ and the density is equal to 
2 



(2 



, , 2 

\v\ < -K 
11-3 



-K < V < 2k 



(4.7) 



Remark 12 The invariant measure for the baker transformation Sp of Example^ is the product 
of the distribution 0/ (1 — A)Ca '^iT'd the normalized Lebesgue measure. Thus, in this example the 
limiting measure for Vn is the distribution 0/ (1 — A)Ca; which may be either singular or absolutely 
continuous. 
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4.2 Graphical illustration of the velocity density evolution with dyadic map 
perturbations 

What is the probabihty density function of /in (Co) defined by Equation 14.21 when is distributed 
according to vl For many maps, including the dyadic map example being considered here, this can 
be calculated analytically which is the subject of this section. 

Let Y be an interval and let v have a density g^. with respect to Lebesque measure. Then the 
distribution of hn{^o) is given by 

Pr{/in(eo) G ^} = Pr{eo G K\A)} = [ g,{y)dy, A G ^(R). 

Jh-\A) 

To obtain the density of /in (Co) with respect to the Lebesgue measure, one has to write the last 
integral as f^gn{x)dx for some nonnegative function gn- If the map /in : 1^ — > R is nonsingular 
with respect to Lebesgue measure, then the Frobenius- Perron operator P^^ : L^([— 1, 1]) — > L^(R) 
for /in exists and gn = Ph„g*- 

Let Y = [—1,1] and let T be the dyadic map. Remember that Pt has the uniform invariant 
density 5*(y) = ^l[_i^i](y). Since hn is a linear function on each interval (^, ^fe^] with constant 
derivative, say h'n, we have 

^ 2"-l 

9n{v) = Ph„g*{v) = —r Yl ^^-^^ 

The derivative satisfies the recurrence equation h'n = A""^ + 2h'^_i, n > 1 and /iq = and is 

equal to ^ — for each n. 

^ A-2 

In Figure 14.11 we show the evolution of the velocity densities gn when T is the dyadic map for 
two different values of A. For A = ^ (left hand panels) the density rapidly (by n = 8) approaches 
the analytic form given in Equation 14.71 On the right hand side, for A = 0.8 the velocity densities 
have, by n = 8, approached a Gaussian-like form but supported on a finite interval. In both cases 
the support of the limiting densities is in agreement with Proposition IHI Figure E21 shows gsiv) for 
six different values of A. 

4.3 r-dyadic map 

Let r > 2 be an integer. Consider the r-dyadic transformation on the interval [0, 1] 

T{y) = ry (mod 1), y G [0, 1]. 

The proof of PropositionlSJcarries over to this transformation when (a^.) is a sequence of independent 
random variables taking values in {0, 1, ... ,r — 1} with equal probabilities, i.e. Pr(afc = i) = ^, 
z = 0, 1, . . . , r — 1, and C is a random variable uniformly distributed on [0,1] and independent of 
the sequence of random variables (a^). Then the limiting measure for Vn is the distribution of the 
random variable 

\ k=0 / 
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Figure 4.1: This figure shows the successive evolution of the densities Qn, n = 2, . . . , 8 of the velocity 
under perturbations from the dyadic map. In the left hand series of densities, A = ^, while on the 

right A = 0.8. The densities when A = ^ rapidly approach the limiting analytic form given 
in Equation 14.71 In both cases, k = 1. 
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A = 0.3 



A = 0.4 



A = 0.55 




Figure 4.2: This figure illustrates the form of the density gs{v) for perturbations coming from the 
dyadic map, as computed from Equation 14.81 and various values of A as indicated, with k = 1. In 
every case the initial velocity density was the uniform invariant density of the dyadic map. 



5 Gaussian Behaviour in the Umit 



r 



0. 



In a series of papers iBeck a,nd Roenstorfll (Il987h . iReckl dlAQOhh . iReckl dlAQOah . iBeckl and 
Hilgers and Bedd ( 19991 ) , motivated by questions related to alternative interpretations of Brownian 



motion, have numerically examined the dynamic character of the iterates of dynamical systems of 
the form 



Un+l = An„ + ^/Tyn 
Vn+l = T{yn), 



A 



--yr 



(5.1) 
(5.2) 



in which T is a 'chaotic' mapping and r is a small temporal scaling parameter. They refer to these 
systems as linear Langevin systems, and point out that they arise from the following 



u 



--fu + y/r'^yi-idit - it). 



Integrating this equation one obtains Equation 15 . 1 1 with u„ = u{nT). 

For situations in which th e map T is selected from t he class of Chebyshev maps [c.f. Equation l2.8l 
and lAdler and RivlinI (|l964h ]. iHilgers and Beckl (l2nOlh have provided abundant numerical evidence 
that the density of the distribution of the sequenc e of iterates itt,;|J^n, for N quite large, may be 
approximately normal, or Gaussian, as A — > 1, and Hilgers and Beckl (|l999^ have provided some of 
the same type of numerical evidence for perturbations coming from the dyadic and hat maps. Our 
results provide the analytic basis for these observations. 

In this section we consider and answer the question when one can obtain Gaussian processes 
by studying appropriate scaling limits of the velocity and position variables. We first recall what 
is meant by a Gaussian process. 
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An R- valued stochastic process {(^{t);t G (0, cxo)} is called Gaussian if, for every integer I > 1 
and real numbers < ti < t2 < ■ ■ ■ < ti < oo the random vector {((ti),... ,C(iO) ^ joint 
normal distribution or equivalently, for all dj £ H, j = 1, . . . , k, the random variable X^j=i djC,{tj) 
is Gaussian. The finite dimensional distributions of a Gaussian process are completely determined 
by its first moment m{t) = E{C,{t)) and its covariance function 

Ki^{t, s) = E{C{t) - m{t)){C{s) - m{s)), s,t > 0. 

If m(t) = 0, t > we say that is a zero-mean Gaussian process. The initial random variable 
C(0) can be either identically equal to zero or can be any other random variable independent of the 
process {C{t);t £ (0,00)}. 

Now we recall the Ornstein-Uhlenbeck theory of Brownian motion for a free particle. The 
Ornstein-Uhlenbeck velocity process is a solution of the stochastic differential equation 

dV{t) = -jV{t)dt + aodw{t), 

where w is a standard Wiener process, and the solution of this equation is 

V{t) = e^^V(O) + ao f e~^^^-'Uw{s). 
Jo 

In other words, V is an Ornstein-Uhlenbeck velocity process if C defined by C{t) = V{t) — e~'^^V{0), 
t > 0, is a zero-mean Gaussian process with covariance function 

2 

Kc{t,s) = ^(e27mm(M) _ i)e-T(*+^). (5.3) 
27 

2 

If the initial random variable V{0) has a normal distribution with mean zero and variance then 
V itself is a stationary, zero-mean Gaussian process with covariance function 

Ky{t,s) = ^e-^\'-< 
27 

Let X{t) denote the position of a Brownian particle at time t. Then 

X{t) = X(0) + / V{s)ds. 
Jo 

In other words, X is an Ornstein-Uhlenbeck position process if rj defined by r]{t) = X{t) — X{0) — 
^— y^y(O) is a zero-mean Gaussian process with covariance function 

2 

K^(t, s) = ^ (27 m.m{t, s)-2 + 2e"^* + 26'^' - e"^!*"'^! - e-^(*+'^)) . 

2 

In particular the variance of T](t) is equal to -^{2^1 — 3 -|- 46""^* — e"^"^*). 

Let h S L'^{Y,l3,i') be such that f h{y)v{dy) = 0. Assume that xq, vq, and ^0 are independent 
random variables on iY,B,u) and ^0 is distributed according to u. The solution of Equation 13.171 
is of the form 

v{t) = e-^(*-"^)i;(nT), t G [nr, (n + l)r), n > 0. 
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We indicate the dependence of x{t) and v{t) on r by writing Xr{t) and Vr{t) respectively. Let 
n = [|:], where the notation [•] indicates the integer value of the argument, for t G [nr, {n + l)r 
substitute A = e~^'^ , and use Equation 13.261 to obtain 



v,{t) = e-^vo + Ke--" e^^'HTHCo)), * > 0, (5.4) 

j=0 

and Eouation 13.301 to obtain 

xr{t) = xo + ^^vo + - V HT^m - e-^* V e^^^hiT^iCo)) ■ (5.5) 

^ U J 

Observe that the first moment of Vr{t) is equal to 

Vr{t)u{dy) = e~^^ / VQ{y)u{dy), 



since the random variables h(T^ {^q)) have a first moment equal to 0. Assume for simplicity that 
J h{y)h{T^ {y))v{dy) = for j > 1 and set = J h?'{y)i'[dy). Since the random variables vq and 
/i(T-^(^o)) are independent, the second moment of VT-{t) takes the form 

jv^tfuidy) = e-'^' [vl{yHdy) + K'e-^^'Y.'''^'' [ ^'^'^'(yX^y) 

= e-'^' J vUyHdy) + a\\-'^'^—^-^^. 
If a and 7 do not depend on r, we have 

Hence the limit of the second moment of Vr{t) as r — > is finite and positive if and only if k 
depends on r in such a way that 

lim — 

r^O r 

is f inite and positive- 
Beck and Roepstorfj (1987) take Kr = y/r from the outset, and claim that in the limit r — > 



the process Vr{t) converges to the Ornstein-Uhlenbeck velocity process when the sequence (hoT^) 
has a so called (/)-mixing property^ on the probability space {Y,B,h'). In fact, the following result 
can be proved. 



sequence of random variables {^j j > 0} is called (^i-mixing if 

JPrM nB) -PrM)Pr(B)| , . , , „ 

hm sup{J ^ ' ^ ' ^ ■.AeT^,be T^+r,, fc > 1} = 

where Ti and J-^„ denote the cr-algebra generated by the random variables £,1, ■ ■ ■ ,(,k and (,k+n, £,k+n+i, ■ ■ ■ respec- 
tively. 
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Theorem 9 Let (y, v) he a normalized measure space and T -.Y he ergodic with respect to 
V. Let h G L'^{Y,B,v) he such that X^^q I I^T i/^l b < and let 

(oo \ 1/2 

j h{yfu{dy)+2Y, I h{y)h{T^{y)Mdy)\ 



n=l 

he positive. Assume that 7 > and 

lim ^ = 1. 

T-»o r 

Then for each vq the finite dimensional distrihutions of the velocity process Vr given hy Equation \5.4\ 
converge weakly as t ^ to the finite dimensional distrihutions of the Ornstein- Uhlenheck velocity 
process V for which V{0) = vq and (Tq = a. 

Proof. By Theorem El we have VT,u{h) = where h = h + f — foT and / = X^^i 'Pj- u^- For 
t > and T > define 

[^] [^] 

j=0 j=0 



Then 

Observe that 



Vr{t)=e'~""vo + Cr{t)+CT{t). 



Cr{t) = Kre-^\e^^^-^ f o Tt^l+i - /) + Kre-^\e-^^ - 1) e^^''/ ° 

j=l 

Hence 

||C.(t)||2<2|K.|e-T*(e^-[^] + 1)11/112, 

and consequently 

llCrWIb <4|k,|||/||2, t>0, r>0. 

This and Lemma 13 imply that the finite dimensional distributions of VT-{t) — e~^^VQ converge weakly 
to the corresponding finite dimensional distributions of a zero mean Gaussian process C with C(0) = 
and the covariance function K(^{t,s) given by Equation 15.31 where cJq = ||/i||25 which completes 
the proof. 

For the corresponding position process we have the following. 

Theorem 10 Under the assumptions of Theorem\^ let V{0) = Vq. Then for each xq the finite 
dimensional distrihutions of the position process Xr given hy Eauation \5.5\ converge weakly as t ^ 
to the finite dimensional distrihutions of the Ornstein- Uhlenheck position process X for which 
X{0)=xo. 

Proof. This follows from Lemma El similarly as the preceding theorem follows from Lemma 

Example 14 Let us apply Theorem|2lto a transformation T : [—1, 1] [—1, 1] and h{y) = y. We 
have VT,uh = when T is the hat map, Eouation 11.41 Then cr^ = 1/3. Thus all assumptions of 
Theorem I2I are satisfied. When T is one of the Chebyshev maps (|2.8|) Sisi we also have VT,uh = 
by Example ITT] and cr^ = 1/2. For the dyadic map (|2.1U|) . the series Y1'^=i'Ptu^ absolutely 
convergent in 1, 1], B{[—1, 1]), u) and is equal to h. This implies that a = \ \2h — hoT\ I2, thus 

in this case cr = 1, which can be easily calculated. Thus all of the numerical examples and studies 
of Beck and co-workers cited above are covered by this example. 
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6 Discussion 



In this paper we were motivated by the strong statistical properties of discrete dynamical systems 
to consider when Brownian motion like behaviour could emerge in a simple toy system. To do 
this, we have reviewed and significantly extended a class of central limit theorems for discrete time 
maps. These new results, presented primarily in Sections 12.31 and 12.41 were then applied in Section 
121 to the Langevin-like equations 

in which the underlying noise r]{t) need not be a Gaussian noise but may be substituted by 



riT) 

n=0 



with a highly irregular deterministic function i^(nr). When the variables /i(^(nr)) are uncorrelated 
Gaussian distributed (thus in fact independent) random variables then the limiting distribution of 
v{nT) is Gaussian. This is need not be the case for the deterministic noise produced by perturbations 
derived from highly chaotic semi-dynamical systems. However (Sectional), in the limit r — > both 
types of noise produce the same stochastic process in this limit, the Ornstein-Uhlenbeck process. 
Finally, in Section |1] we have illustrated all of our results of Section |21 for the specific case of 
perturbations derived from the exact dyadic map. 

The significance of these considerations is rather broad. It is the norm in experimental observa- 
tions that any experimental variable that is recorded will be "contaminated" by "noise" . Sometime 
the distribution of this noise is approximately Gaussian, sometimes not. The considerations here 
illustrate quite specifically that the origins of the noise observed experimentally need not be due 
to the operation of a random process (random in the sense that there is no underlying physical law 
allowing one to predict exactly the future of the process based on the past). Rather, the results 
we present strongly suggest, as an alternative, that the fluctuations observed experimentally might 
well be the signature of an underlying deterministically chaotic process. 



Acknowledgments 

This work was supported by MITACS (Canada) and the Natural Sciences and Engineering Research 
Council (NSERC grant OGP-0036920, Canada). The first draft of this paper was written while 
the second author was visiting McGill University, whose hospitality and support are gratefully 
acknowledged. We thank Radoslaw Kamihski for providing the computations on which Figure ESI 
was based. 



A Appendix: Limit Theorems for Dependent Random Variables 

This Appendix reviews known general central limit theorems from probability theory which we 
can use directly in the context of noninvertible dynamical systems. With their help we then prove 
several results which we have used in the preceding Sections. 
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Consider random variables arranged in a double array 

Ci,i) Ci,2i • • • 1 Ci,fci 

(2,1,(2,2, ■■■ ,C2,k2 (1-1) 
Cn.l) Cn,2, ■ ■ ■ , Cn,kn 



with A;„ — > 00 as n — > CO. We shall give conditions which imply that the row sums converge in 
distribution to a Gaussian random variable aN{0, 1), that is 

Y,Cn,^^''aNiO,l). (1.2) 
1=1 

We shall require the Lindeberg condition 

i™oE ^(Cn,il{|C„.|>4) = for every e > 0. (1.3) 

^ j=i 

If independence in each row is allowed, then we have the classical Lindeberg- Feller theorem. 

Theorem 11 Theorem 7.2.1)) Let the random variables Cn,i, ■ ■ ■ ,Cn,kn ^6 indepen- 

dent for each n. Assume that E{(^n.i) = and X^f=i -E'(Cn i) = ^, n > 1. Then the Lindeberg 
condition holds if and only if 

Y.Cn,^ -'iV(0,l), 
i=l 

and 

for all (5 > 0, max Pr{|C„,i| > 5} ^ {). (1.4) 

l<i<kn ' 

If one considers a map T on a probability space {Y, B, u) which preserves the measure and defines 
Cn,i to be o for i = 1, . . . , n with / measurable, then Condition II .41 holds. This is because 

Pr{|C„,.| >5} = u{{y G Y : \f{T'-\y))\ > 6^^}), 

and {y £Y : \ f{T'^\y))\ > = T^'+Hiy e ^ : l/(y)l > ^V^})- Thus by the invariance of 

this leads to 

max Pr{|Cn,i| > 6} = u{{y G Y : \f{y)\ > 6^}) ^ 0. 
i<i<fe„ 

Similarly, if one takes a square integrable /, then the Lindeberg condition 11.31 is satisfied. Indeed, 
we have 

1 ' e2rrri-ii 



^(lCn,.ri{|C„,|>e}) = -=2 / f\T^-\y)Hdy) 



and by the change of variables applied to T* ^ this reduces to 

- [ f\y)v{dy). 
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Hence 

n ~ 

^£(|Cn,pl{|C„„|>e})= / fivMdy), 
~i J{\f\>V^^} 

which converges to by the Dominated Convergence Theorem under our assumption that is 
integrable. 

Since our random variables are dependent, we cannot apply the above theorem. Instead, we 
use the notion of martingale differences for which there is a natural generalization of Theorem 1111 
Moreover, additional assumptions are needed, as one can easily check that if T is the identity map, 
then 

V" , V ™ 

1=1 

which can not be convergent to a Gaussian random variable. 

We first recall the definition of conditional expectation. Let a probability space (0, T, Pr) be 
given and let ^ be a sub-cr-algebra of J-. For C G L^(r2, J^,Pr) there exists a random variable 
^'(Cl^)) called the conditional expected value of C, given Q, having the following properties: it is Q 
measurable, integrable and satisfies the equation 

/ E{C\G){u:)Vi{duj) = [ C{^)FiidLo), A eg. 

J A J A 

The existence and uniqueness of E{(^\g) for a given ^ follows from the Radon-Nikodym theorem. The 
transformation ^ i— > E{(^\g) is a linear operator between the spaces L^{i^, J-,Pv) and L^(0,^,Pr), 
so sometimes it is called an operator of conditional expectation. 

Let {Cn,i : 1 < ^ < knjTi- > 1} be a family of random variables defined on a probability space 
(r2,.F, Pr). For each n > 1, let a family {J^n.,j : i > 0} of sub-cr-algebras of J- be given. Consider 
the following set of conditions 



(i) E{Cn,i) = and E{Cli) < oo, 

(ii) J' n,i~l ^ ri,i, 

(iii) Cn,i is Tn,i measurable. 



(1.5) 



(iv) E{Cn,i\Fn,i-i) = for each I <i < kn, n>l. 



A family {J-n,i,Cn,i : 1 < ^ < kn,n > 1} satisfying conditions (i)-(iv) is called a (square integrable) 
martingale differences array. 

The next theorem is from iBilhngslevI (Il99,5h . 



Theorem 12 \Billinaslei\ h99A . Theorem 35.12)) Let {Cn.i : 1 < ^ < kn,n > 1} he a martingale 
difference array satisfying the Lindeberg condition \1.3l If 



where a is a nonnegative constant, then 

Y,Cn,i^'' oN{^,l). (1.6) 
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If the limit in Condition ll.Sl is a rand om variabl e inste ad of the constant o"^, we obtain conver- 
gence to mixtures of normal distributions lEag lesoni dl 97,4 Corollary p. 561). 

Theorem 13 Let {Cn,« : 1 < ^ < kn,n > 1} be a martingale difference array satisfying the Linde- 
berg condition M.,"^ If there exists an J^^o = {^^^iJ^nfl-'measurable, a.s. positive and finite random 
variable rj such that 

Y,E{Ch\^n,-i)^''v, (1-7) 

i=l 

then X]f=i Cn,i is convergent in distribution to a measure whose characteristic function is (p{r) = 
E{exp{-^r'^ri)). 

The above result shows that to obtain a normal distribution in the limit a specific normalization 
is needed and we have the following 

Theorem 14 \Gaenssler and Jooi \l99A . Theorem 3.6)) Let {Cn.i : 1 < i < kn-,fi > 1} be a 
martingale difference array and let r] be a real-valued random variable such that Pr(0 < r] < oo) = 1. 
Suppose that 

lim E{ max |Cn,i|) = (1.8) 

n— >oo l<i<k„ 

and 

ECn.-^^- (1-9) 

i=l 

If T] is J^n^i-measurable for each n and for each 1 < i < kn, then 



Remark 13 If r] in Condition M.IA is constant, rj = a"^ , then the conclusion of Theorem is 
equivalent to 

Y,Cn,^ aiV(0,l). 

i=l 

Note also that Condition M.t^ is implied by the Lindeberg condition. 

The next result gives conditions for moment convergence in Theorem 1121 

Theorem 15 \Hal\ Iil974. Theorem MTeicheA 1198 A . Theorem 3)) Let {Cn,i : 1 < « < A;„, n > 1} be 
a martingale difference array with ^^=1 ^iCni) ~ where o" > 0. Suppose that for p > 1 

J^E|Cn,.P^-0 and E\Y,E{Cli\^n,-i)-o^\P^^. (1.10) 

i=l 1=1 

Then 

kn 

lim E\y2CnA^P = E\NiO,a^)\^P. 

i=l 
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Let {Y, B, u) be a normalized measure space and T -.Y he a, measurable map such that T 
preserves the mcasurG v. Rcca-ll from S6ction 01 th.6 rclcttion 12 . 3l between tlie transfer operator I^t^v-^ 
the Koopman operator and the operator of conditional expectation which gives 

VT,uoUTf = f and UToVT,uf = E{f\T~HB)), f € lHy,B,u). 

Lemma 1 Let (Y, B, u) he a normalized measure space, and T : Y ^ Y be a measurable map such 
that T preserves the measure u. Let {c„^j : 1 < "i < kn,n > 1} be a family of real numbers and 
h G L'^{Y^B,v). Suppose that VT,uh = 0. Then 

Cn,i — Cn^i h oT " ^ 1^^^!^ ^ru Cn,i — 0; ^ ^ ^rxj 

with 

Tn,i = T~^'^+\B), 0<i<kn, and T^a = B, i > K, n > 1 
is a martingale difference array and ifcni = ^ ^ i ^ kn then the following hold 

(i) Lindeberg condition ] 1. 'A 

(ii) Conditions and Al.fA with rj = E{h^\T) where I is the a-algebra of all T -invariant sets; 

(iii) Condition \1.5\ with o"^ = J K^dv provided that T is ergodic; 

(iv) Condition M.im for every p > 1 provided that T is ergodic and h € L°"{Y,B,i'). 

Proof. To check conditions (ii), and (iii) of the definition of a martingale difference array, observe 
that T~^~^{B) C T~^{B) and h o is T~^{B) measurable. The Koopman and transfer operators 
for the iterated map are just the j*'^ iterates of the operators Ut and Vt,v From this and 
Equation 12 .HI we have Vj^JJ^h = h and 

E{h o T^\T~^~^{B)) = C/^+V^+\/i o T^) = U^-^^VT,^h. 

Since VT,uh = 0, we see that E{h o T^\T~^~^{13)) = for j > which proves condition (iv). 
The Lindeberg condition reduces, through a change of variables, to 



k 

~l J{\h\>VKe} 



but /i^ is integrable and the Lindeberg condition follows. To obtain Condition 11.71 use Equation 
12.31 change the order of summation 

Y,E{C,\:Fn,-i) = -^£;(/i2or'="-|r-^"+-i(fi)) 

1=1 1=1 

^'^ i=i 

1=1 i=0 

and apply Birkhoff's ergodic theorem to the integrable function UT'PT,uih'^) to conclude that this 
sequence is convergent to E{UT'PT,v{h'^)\1) almost everywhere (with respect to v), and consequently 
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in probability. Since UtVtA^'^) = E{h?\T~^{B)) and I C T~^{B), we have E{UtVtA^'^)\^) = 
E{K^\I). Similarly, Condition 11.91 follows from the Birkhof ergodic theorem. In addition, if T is 
ergodic, then rj is constant a.e. and is equal to J h?di'. Since h? G LP{Y,B, v) and p > 1, we have 

i=i 

and UT'PT,u{h^) £ LP{Y,B,i'). By the ergodic theorem in spaces we get 

hji -i — ■ 1 p 

1=1 1=0 

in LP{Y,B,i'), but J UT'PT,u{h'^){y)v{dy) = f h?{y)i'{dy) and the proof is complete. 

We now turn to the FCLT for {h o T*)j>o. Let (YjB,!^) be a normalized measure space and 
T : y — > y be ergodic with respect to u. Let h £ L'^{Y,B,i') and a = \\h\\2 > 0. We define a 
random function 

[nt]~l 

i^nit) = ^ V /i o r for t G [0, 1] 

(where the sum from to —1 is set to be 0). Note that ijjn is a right continuous step function, a 
random variable of -D[0, 1] and V'n(O) = 0. 

Lemma 2 If VT,uh = 0, then the finite dimensional distributions of ipn converge to those of the 
Wiener process w. 

Proof. To show that the finite dimensional distributions of ipn converge to the corresponding finite 
dimensional distributions of w we use the Cramer- Wold technique. If the ci, . . . ,Cfc are arbitrary 
real numbers and io = < ti < . . . < tfc < 1, we put 



Cn.i 



'^h o T""*, n - [ntj] <i<n - [ntj^i], j = l,...,k 



0, 1 < i < n — [ntk] and tfc < 1 

, 0, i > n 



Observe that 

n k 



i=l j=l 

By Lemma nCn,i is a martingale differences array and we will verify the conditions of Theorem 1121 
For the Lindeberg condition note that 

k ^2 

i=l j=l 

and as a finite sum of sequences converging to it is convergent to 0. For Condition 11.51 observe 
that 

n k 2 k 

EEicl\T-^'~HB)) = j2j^ E umh')^^^c]it,-t,^,) 

i=l j=l i=[ntj_i]+l j=l 
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by the ergodicity of T, and the fact that = J K^du. Therefore, by Theorem 1121 

n 



i=l 



^ 5]c2(t,-t,_i)iV(0,l). 



Thus X]j=i Cj('0n(*j) — V'n(^i-i)) converges to the Gaussian distributed random variable with mean 
and variance X]j=i'^^(*i ~ ^i-i)? t>ut this is the distribution of l^j=i Cj (u)(tj) — tt;(tj_i)) which 
completes the proof. 

Lemma 3 IfVT,uh = 0, then Condition \2.11A holds for each positive e. 
Proof. Define a martingale difference array 

f 1 /loT"-*, l<i<n 

/- . — ) ay/n ' — — 

I U, z > n, n > 1. 

Let also Cn,o = and = J2i^o Cn,i, t G [0, 1], n > 1. We have 

We first observe that 

sup \lpn{s) - ^pn{t)\ < sup | V'n (•?) " ^n(0 1 
|t-s|<(5 |t-s|<'5 

< 4 sup sup - Vn(A;(5)|. 

k kS<t<(k+l)6 

This gives 

u{ sup IV'n('S) - -^nC*)! > e) < SUp jV'nW " ^n(H)| > 7)- 

\t-s\<S 1^^^^ kS<t<{k+l)5 4 

Now applying Lemma |21 and arguments similar to those of Brownl (1971, pp. 64-65), one can 
complete the proof. 

For the next results we need the following 

Lemma 4 Let {zi)i>i be a sequence of real numbers such that 

lim ^^1^ = z. 

n— >oo n 

Then 

n— >oo sr^'^n i 

for every sequence of integers kn > 1 and every sequence of real numbers a„ satisfying 



lim kn = oo, lim a„ = 1, and lim a^" ^ 1 (l-H) 

n— >oo n— ♦oo n— >oo 

and either a„ > 1 or < < 1 /or all n > 1. 
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Proof. The Abel method of summation, 



m—l 



^ Cibi = Cm'^bi - - q) ^ bi 



i=k 



i=k i=k 



j=k 



can be used to write 



kfi 1 



^^n = - 5Z ^^'^'^ ~ Zj - iz) + 0^"(^ Zi - K.i 



i=l i=l 
Fix e > and let no be such that 



i=l 



Z^i=i Zi 



m 



< e for m > uq. 



i=l 



Suppose that a„ > 1 for all n > 1. The other case is proved analogously. Combining these yields, 
for kn> riQ, 



i=l "n 



< 



1=1 "n 



1 = 1 



Letting n ^ oo we see that the first term on the right goes to zero, while the second term goes to 
e times a constant not depending on e, which completes the proof. 

Theorem 16 Let (y, v) be a normalized measure space and T : Y ^ Y be ergodic with respect 
to V. Let {kn), (cin) be sequences satisfying Condition \l.ll\ Let c G R and let (cn) be a sequence 
of real numbers such that 

lim knC^ = (?. 

n— too 

If h £ L'^{Y, B, u) is such that T'T,uh = 0, then 



CnY,a}^hoT^^ cTiV(0,l), 



i=l 



where a 



6\a'-l) 
In 



2 an- 



d a = lim aj" . 



Proof. From Lemma ^ it follows that Cn,i = Cna^~^^ */i o T^^+i « jg ^ martingale difference array 
and that 



lim 



1 



We shall apply Theorem I14L We have 



E{ max \hoT^\) = 0. 



max |Cn,i| < |cn| niax(a^", 1) max |/ioT*|, 

l<i<kn ' l<i<kn 



so 



£'( max |Cn,j|) < V ^nlcnl max(a„", = max {hoT]). 

l<i<k„ y kr, 0<i<kn 
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Letting n — > oo we see that Condition 11.81 holds. To verify Condition II .91 note that 

i=l i=l i=l l^i=l 

Therefore Birkhoff's ergodic theorem, Lemma HI and the fact that 

hm <2^a^* = -\ ^ 

n— >oo ma^ 

i=l 

complete the proof. 

Corollary 5 Under the assumptions of Theorem Mb\ if for h G L/^iY^B^v) the series Yl'n=o'^T 
is convergent in L'^(Y,B,i'), then 



CnJ^aihor^'' aiV(0,l), 



"'^ere a = \/ ^^^^ \\h + f - f o T\\2 and f = Zn=i ^l^h 



Proof. Theorem IHl implies VT,u{h + f — f °T) =0. Thus 

cn <ih + / - / o T) o r ^'^ aNiO, 1) 

i=l 

by Theorem 1161 Therefore it remains to prove that 



cnJ^<(/or-/)or ^^0. (1.12) 

i=l 



Observe that the left-hand side of Equation 11.121 is equal to 

cn{ai-f o T^-+^ -foT) + cn{a-' - 1) aj,/ o T\ 



i=l 



Since — > as n — > oo, the first term converges in probability to 0. From Lemma^Jand Birkhoff's 
ergodic theorem it follows that 



i=l "n 



f{y)v{dy). 



Therefore the sequence 

c.(a-i - 1) 5^ a^J o r = c.(l - a',r) ^^^C , 

i=l 1^1=1 

is also convergent in probability to 0, which completes the proof. 



59 



Remark 14 Note that we can conclude from Theorem \16\ that 

i=m 

where m > is any fixed integer, because Cn goes to zero and a„ to 1 as n ^ oo, so the difference 

CnJ2<^°T'-Cn^al,hoT\ 

1=1 i=m 

which is either equal to Cnh or Cn ^^Li ° converges in probability to zero. 

Lemma 5 Let {Y,B.,u) be a normalized measure space, T : Y ^ Y be ergodic with respect to v, 
and J ^ be a constant. Let k,-, r > 0, be such that 

lim ^ = 1. 

r^O T 

Lf h £ L'^(Y,B, v) is such that T'T,uh = 0, then the finite dimensional distributions of the process 
defined by 

[^] 

converge weakly as t ^ to the corresponding finite dimensional distributions of the zero-mean 
Gaussian process C for which ^(0) = and 

EC(t)C(s) = Ml(e27min{t,s) _ l^^-jit+s) ^ t,S>0. 
27 

Proof. To prove the convergence of the finite dimensional distributions of Ct to the corresponding 
finite dimensional distributions of the Gaussian process it is enough to prove that for any / > 1, 
real numbers < ti < ... < t; < 00 and di, ...,di the distribution of X]j=i djC{tj) is Gaussian and 

that X]j=i djCr{tj) converges in distribution as r ^ to Yl'^j=i djC{tj). 

We consider first the case of / = 1. It follows from Theorem 1161 that for t > the distribution 
of Cr{t) converges weakly as r — > to a Gaussian random variable. To see this let r„ be a sequence 
going to zero as n — > 00. Take A;„ = [— ], On = e'^'^", Cn = Kr^e""^*, and observe that 



lim kn = 00, lim = 1, lim a^" = e"^* and lim fe^c^ = te 

n— »oo n— >oo n—KX) n— >oo 

and — — iJ'^i-yt — — = ^~2^ ^ . The theorem then implies that 



■27* 



where i7q = J h'^{y)v{dy) and ^0 is distributed according to v. Consequently Cr(0 C(0 as r ^ 0, 
where C,{t) is a Gaussian distributed random variable with mean and variance given by 



27 



i(l_e-27i)^ t>0. 
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Note that Ct(0) = Hrh- Since limT-_+o = 0, we also have Ct(0) ^ as r ^ 0. 

We next consider the case of / = 2. The case of arbitrary / is deduced analogously from Theorem 
[TH Let ti < t2 and di, d2 be given. Let r„ be a sequence going to zero as n — > oo. Set kn,i = [^], 
kn,2 = kn = kn,2 + Ij and observe that < fcn,2 for all n sufficiently large. Define 

0, otherwise. 



Then we have 



Observe that 



fc„ kn,2 

j=0 

As in the proof of Theorem 1161 we check that Theorem 1141 applies to {r]n,j ■ ^ ^ j ^ kn^n > 1} and 
conclude that 

C?lCr„(tl)+d2Cr„(t2) ^'^ CTiV(0, 1) 

where = ^(^1(1 - e-^T*^) + 2(i2die-^(*2+ti)(e27ii _ i) + ^2^^ _ £-27*1)). Since fjiV(0, 1) is 

the distribution of (iiC(ti) + d2C(i2) and ^(C(ti)C(i2)) = Mle-7(t2+ti)(e27ti _ 1)^ the proof of the 
lemma is complete. 

Lemma 6 Let {Y,B^v) he a normalized measure space, T : Y ^ Y he ergodic with respect to v, 
and J ^ be a constant. Let Kr, t > 0, he such that 

K-2 

hm ^ = 1. 

T-»o r 

Lf h ^ L'^{Y,B,u) is such that VT,vh = 0, then the finite dimensional distributions of the process rjr 
defined hy 



7j^{t) = -L ^(1 - e^^^-'-*))/! o T-J', t > 0, T > 



converge weakly as t ^ to the corresponding finite dimensional distributions of the zero-mean 
Gaussian process r] for which r]{0) = and 

Er]{t)rj{s) = ^^(27 min(t, s) - 2 + 26"^* + 2e-^' - e-Tl*"^! - e-^(*+*)) 

for t,s > 0. 

The lemma follows from Theorem 1141 in a similar fashion as the preceding lemma. 
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